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1 . 0  INTRODUCTION 

The  purpose  of  this  study  is  to  document  the  sensitivity  of 
multivariate  statistical  interpolation  to  the  values  of 
statistical  parameters,  misspecif ication  of  these  parameters,  and 
to  missing  observations.  It  is  hoped  that  this  information  will 
enable  practitioners  to  concentrate  on  the  appropriate 
specification  of  crucial  parameters  while  avoiding  agony  over 
those  that  have  relatively  small  effect.  The  results  presented 
here  point  out  some  interesting  dependencies  which  I  believe  have 
not  been  previously  published.  Prior  studies  upon  which  this 
work  builds  are  Franke  (198b  and  1988)  and  especially  Franke,  et 
al.  (1988).  It  was  also  influenced  by  Seaman  (1983).  These 
previous  studies,  however,  involved  only  univariate  objective 
analysis  schemes.  While  many  of  the  results  can  be  expected  to 
carry  over  in  a  similar  way,  the  influence  of  wind  observations 
on  the  expected  error  over  various  observation  sets  is  useful  and 
interesting.  Further,  I  believe  that  the  graphical  presentation 
of  the  results  used  here  makes  it  easy  to  discern  the  important 
parameters . 

A  brief  overview  of  multivariate  statistical  interpolation 
and  the  method  used  to  calculate  expected  errors  is  given  in 
Section  2.  The  details  of  the  necessary  calculations  are  not 
given  explicitly  there,  but  are  available  in  Appendix  2,  which 
gives  a  listing  of  the  subroutine  used  to  evaluate  the  covariance 
matrices.  Section  3  gives  a  discussion  of  the  various  parameters 
and  the  values  over  which  they  were  varied  during  the  study. 
Section  4  contains  an  analysis  of  the  results.  Because  of  the 
plethora  of  data  which  was  generated  it  is  difficult  to 
comprehend  the  important  details  in  tabular  format,  and  therefore 
the  Information  is  incorporated  into  a  few  graphs  which  enable 
one  to  easily  ascertain  the  sensitivities  of  the  scheme  to  the 
various  parameters.  Tables  are  included  in  Appendix  1  for 
completeness,  but  it  is  expected  that  few  readers  will  find  them 
necessary.  Section  5  contains  the  concluding  remarks. 
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2.0  MULTIVARIATE  STATISTICAL  INTERPOLATION 

Statistical  Interpolation  (SI)  is  in  use  at  many  of  the 
world's  Numerical  Weather  Prediction  (NWP)  facilities,  including 
the  U.S.  Navy's  Fleet  Numerical  Oceanography  Center  (PNOC).  The 
scheme  has  its  roots  in  the  work  of  Kolmogorov  and  Weiner,  and 
was  first  developed  for  meteorological  applications  by  Gandin 
(1965).  More  recently  it  has  been  applied  in  multivariate  form 
by  Schlatter  (1975),  Schlatter,  et  al.  (1976),  Bergman  (1979), 
and  Lorenc  (1981). 

In  its  perfect  form  (all  parameters  known),  SI  delivers 
estimates  of  a  field  with  minimum  mean  squared  error  over  a 
certain  ensemble  of  realizations  that  satisfy  normality  and 
stationarity  of  the  underlying  stochastic  process.  For 
convenience,  isotropy  is  usually  assumed,  and  for  meteorological 
problems  a  zero  mean  is  assumed  (although  a  nonzero  mean  can  be 
accounted  for  in  more  than  one  way) . 

In  the  multivariate  formulation,  the  dependent  data  consists 
of  related  variables,  which  in  our  case  we  assume  to  be  the 
errors  in  the  background  pressure  height  and  wind  fields,  H,  U, 
and  V,  as  related  through  the  geostrophic  relationship: 


U  =  k^H^,  V  -  kjHy  . 

Here  the  values  of  kj^  and  are  dependent  on  the  latitude.  The 
SI  equations  are  applied  to  the  background  error,  obtained  by 
forming  the  difference  between  the  background  (normally  the  NWP 
values,  interpolated  to  the  observation  locations)  and  the 
observed  values.  The  assumption  of  normality  implies  the  minimum 
variance  (or  least  mean  squared  error)  predictor  is  a  linear 
combination  of  the  data  values.  Then  construction  of  the  weights 
leads  to  solution  of  a  system  of  equations  whose  coefficient 
matrix  is  the  matrix  of  spatial  covariances  and  cross-covariances 
for  the  variables.  Because  of  the  assumed  relationship  between 
the  variables,  the  wind  error  covariances  and  the  cross¬ 
covariances  are  determined  when  the  height  error  spatial 
covariance  function  is  known.  The  exact  relationship  is  given  in 
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Franke,  et  al.  (1988),  as  well  as  in  other  references  above. 

The  SI  equations  have  been  derived  in  numerous  places  (see 

e.g.,  Schlatter  (1975),  Lorenc  (1981),  ThiSbaux  (1985),  and 

Thi6baux  and  Pedder  (1987))  and  are  repeated  here  only  for 

completeness.  Let  O  =  (o.)  „  represent  N  measured  values, 

with  corresponding  independent  observation  error  variance 
.  2 

(diagonal)  matrix  Z=  {Cj},  B=  {b^}  the  corresponding  vector  of 
background  values,  and  C={c.  „  the  spatial 

covariance  matrix  for  the  background  errors.  Then,  letting 
C,^={c- . } ,  „  denote  the  vector  of  covariances  between  the 

error  in  the  variable  at  the  location  at  which  it  is  to  be 
analyzed  and  the  background  errors,  the  following  equation  holds 
for  the  analyzed  value,  a^: 

a^  =  bg  +  cJ(C+S)''^(0  -  B)  . 

If  the  statistical  parameters  are  known  precisely  (whereupon 
SI  becomes  Optimum  Interpolation,  or  Ol),  then  the  expected  error 
variance  for  the  estimated  variable  at  P^  is  given  in  the  usual 
least  squares  form, 

2  2  T 

‘'o  ‘^b  "  ci(C+E)  , 

2  . 

where  is  the  variance  of  the  background  error.  Considering 
the  case  where  the  statistical  parameters  are  not  known  exactly, 
the  analyzed  value  becomes 

Sq  =  bg  +  cJ(C+S)“^(0  -  B)  , 

where  the  tilde  overbar  signifies  assumed  inexact  values.  Notice 
that  the  equation  is  exactly  the  same  except  that  assumed  values 
replace  the  exact  values  for  covariances.  The  expected  mean 
squared  error  is 

ol  =  -  2cJ(C+£)“^Cq  +  cJ(C+Z)"^(C+E)(C+2)"^Cq  . 
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These  equations  are  given  in  Seaman  (1983),  where  the  vector 
Wq  =  (C+S)  is  interpreted  in  the  usual  statistical 
fashion,  as  the  weight  vector  for  the  observed  minus  background 
values.  This  also  simplifies  the  writing  of  the  error  expression 
and  indicates  a  computationally  efficient  algorithm. 

The  expected  error  will  be  computed  at  a  number  of 
locations,  in  our  case  on  a  7x11  grid  of  points.  The  matrix 
(C+S)  is  symmetric  and  positive  definite,  so  the  computation  of 

W  is  accomplished  by  performing  a  one  time  Cholesky 

^  ~  ~  .  T 

decomposition  of  (C+S)  into  LL  ,  followed  by  forward  and 

backward  substitutions  to  find  W  as  the  solution  of 
T  ~  .  .  .  ^ 

LL  Wq  =  Cq  .  This  IS  an  important  concept  since  the  number  of 
equations  is  the  number  of  observations  (up  to  108  here) ,  and 
varies  over  77  locations  in  our  case.  In  actual  practice 
would  vary  over  several  locations  for  the  same  observation  set 
when  a  block  analysis  scheme  is  used. 


3.0  SETTING  FOR  THE  STUDY 

An  empirical  study  such  as  this  one  requires  a  number  of 
compromises  concerning  the  range  of  parameters  permitted.  In 
previous  studies  (Franke,  1985,  and  Franke,  et  al.,  1988)  a  set 
of  three  grids  and  corresponding  observation  locations  based 
approximately  on  the  radiosonde  network  over  the  United  States 
and  the  Atlantic  Ocean  were  used.  The  study  here  is  also  at  a 
single  level.  One  difference  from  my  previous  studies,  which  has 
a  rather  minor  influence,  is  that  distances  are  measured  in 
meters  (approximately;  I  did  not  use  the  exact  geodesic 
distance)  rather  than  degrees.  The  formula  used  for  the  distance 
(also  used  by  Schlatter  (1975))  between  two  points  (0  = 
longitude,  ip  =  latitude)  is 

dfj  =  r^[(<p.-«p.  )^  +  (0.-e^)^cos^(^l^)]  . 

While  computed  distances  are  (relatively)  quite  different  in  the 
upper  latitudes,  the  overall  influence  on  the  expected  mean 
squared  errors  is  rather  small  over  the  regions  we  consider. 
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Figures  1-3  show  the  2.5°  grids  and  observation  locations 
used  in  this  study.  The  +'s  indicate  the  grid  points,  while  the 
small  squares  and  circles  represent  the  observation  locations. 

The  shaded  circle  represents  the  "missing  observation"  for  tests 
of  sensitivity  of  the  analyzed  values  to  one  missing  observation 
(both  missing  wind  only  and  missing  height  and  wind).  The  open 
circles  represent  additional  missing  observation  locations  for 
tests  of  sensitivity  to  many  (one-half)  missing  observations. 

Each  grid  is  a  7x11  grid  of  points,  taken  to  be  the  interior  grid 
points  of  a  9x13  grid  containing  the  observation  locations.  The 
Middle  United  States  (MUS)  grid  has  36  observations,  the  East 
Coast  (EC)  grid  has  25  observation  locations,  while  Middle 
Atlantic  (MA)  grid  has  only  3  observations. 

The  assumed  correlation  function  for  the  height  errors  is 
the  specialized  second  order  autoregressive  function  (SOAR)  that 
seems  to  be  quite  stable  with  respect  to  variations  in  the 
parameters  while  embodying  enough  parameters  to  fit  historical 
data  (Franks,  et  al.,  1988).  As  noted  by  Thi6baux,  et  al. 

(1986),  Balgovind,  et  al.  (1983),  the  SOAR  seems  to  have  an 
affinity  for  meteorological  data.  The  form  we  used  is 

C(s)  =  (l-A) (l+as)e“^®  +  A, 

where  s  denotes  the  distance,  and  a  and  A  are  parameters  which  in 
practice  are  determined  by  a  fitting  process.  In  the  previous 
section  c.  .=C(d.  .  )a^,  where  d .  .  is  the  distance  between  the  i 
and  j  observation  points.  The  study  conducted  here  made  use  of 
five  different  sets  of  parameter  values,  as  shown  in  Table  1,  and 
also  depicted  in  Figure  4.  The  "nominal”  correlation  function  is 
considered  to  be  #4,  which  approximates  the  Bessel  function  curve 
of  Lonnberg  (1982)  closely,  and  also  corresponds  closely  to  the 
decay  rate  used  by  Lorenc  and  Mammon  (1988).  Varying  from  one 
correlation  function  to  the  next  makes  it  possible  to  determine 
whether  the  additive  constant  or  the  decay  rate  constant  is  the 
more  critical. 
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Cor# 

A 

a 

#m  to  #m+l 

1 

0.0 

5.000x10“® 

2 

0.0 

3.082X10“® 

a 

decreases 

3 

0.15 

3.082X10“® 

A 

increases 

4 

0.2722 

3.082X10“® 

A 

increases 

5 

0.2722 

2.188x10“® 

a 

decreases 

Table  1:  Correlation  Functions 

The  nominal  value  for  the  standard  deviation  of  the 
background  height  error  is  30  m.  The  correlation  function  then 
determines  the  variance  of  the  wind  field  errors,  and  the 
standard  deviations  for  wind  errors  are  given  in  Table  2. 

Because  the  values  of  and  «5®pend  on  the  latitude,  the  rms 
errors  over  a  given  grid  depend  on  the  grid.  All  expected  wind 
errors  were  calculated  relative  to  their  value  at  a  given  point, 
and  the  rms  values  of  these  are  given  in  the  tables  and  the 
figures,  different  than  computing  the  expected  rms  values  over 
the  grid  and  then  comparing  this  with  the  rms  background  error 
over  the  grid. 

The  nominal  values  assumed  for  the  observation  errors  were 
10  m  for  heights  and  l.O  m/sec  for  wind  components.  For  01  these 
values  were  varied  over  the  values  0,  5,  10,  20,  and  »  m  for 
heights,  and  0,  0.5,  1.0,  2.0,  and  «  m/sec  for  the  winds.  The 
expected  rms  error  was  not  computed  for  all  combinations  of  these 
values;  the  <»  values  imply  no  measurement  of  that  variable  and 
computations  were  only  performed  with  the  nominal  value  of  the 
other  measurement  error,  for  example.  These  calculations  were 
perfomned  with  01  to  show  the  effect  of  no  observations  of  a 
particular  variable.  With  no  wind  observations,  the  process 
collapses  to  the  usual  univariate  01  scheme. 

In  order  to  assess  the  effect  of  missing  observations,  four 
"missing  observation"  computations  were  performed.  Foi  one 
observation  point  in  each  grid  it  was  assumed  that  the  wind 
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observations  were  not  available;  then  it  was  assumed  that  the 
height  observation  at  the  same  location  was  not  available.  For 
about  one-half  of  the  observation  points  it  was  assumed  that  the 
wind  observations  were  not  available;  then  it  was  assumed  that 
the  corresponding  height  observations  were  not  available. 

The  computations  of  expected  rms  errors  for  the  statistical 
schemes  did  not  cover  such  a  wide  range  of  assumed  observation 
error  values.  Variations  of  the  assumed  observation  errors  and 
spatial  correlation  function  number  were  only  carried  out  to  the 
adjacent  value  (assumed  observation  error  changing  by  a  factor 
of  .5  or  2,  and  correlation  function  number  only  to  the  adjacent 
number ) . 


Cor# 

grid  1&2 

grid  3 

1 

16.05 

14.52 

2 

9.89 

8.95 

3 

9.12 

8.25 

4 

8.44 

7.64 

5 

5.99 

5.42 

Table  2;  Background  wind  errors,  m/sec 


4.0  ANALYSIS  OF  THE  RESULTS 

The  easy  detection  of  the  sensitivity  to  parameters  on  which 
a  process  depends  is  sometimes  clouded  by  the  mass  of  information 
available.  In  the  following,  I  believe  the  plotting  scheme 
adopted  enables  readers  to  easily  detect  critical  parameters. 

In  this  study  there  are  four  important  parameters  that  are 
varied  for  each  of  the  three  grids;  The  height  observation 
accuracy,  the  wind  observation  accuracy,  the  additive  constant  A 
in  the  correlation  function,  and  the  decay  rate  a  (the  reciprocal 
can  be  thought  of  as  some  measure  of  "correlation  distance")  in 
the  correlation  function.  In  this  study  the  correlation  function 
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paramet'^i.s  have  not  been  varied  except  one  at  a  time,  in  the 
manner  noted  in  the  previous  section. 

We  first  consider  the  sensitivity  of  01  to  the  accuracy  of 
the  observations.  The  results  for  the  error  in  the  analyzed 
height  field  for  13  different  combinations  of  observation 
accuracies  and  5  correlation  functions  are  shown  in  Figures  5-7 
for  the  MUS,  EC,  and  MA  grids,  respectively.  Note  that  only  the 
integer  abscissa  values  have  meaning,  and  the  points  are 
connected  only  to  enable  one  to  more  easily  see  the  effects  of 
changing  parameter  values.  As  one  moves  to  the  right  on  the 
scale,  the  trend  is  toward  less  observation  accuracy,  with  wind 
observation  accuracies  first  decreasing,  while  the  height 
observation  accuracy  more  slowly  decreases.  From  the  plots,  it 
is  seen  that  significant  increases  in  01  height  analysis  errors 
occur  at  abscissae  2,  5,  9,  and  13,  while  between  there  is 
generally  some  increase,  but  much  smaller  in  magnitude.  Table  3 
lists  the  values  of  the  observation  accuracies,  and  it  is  seen 
that  for  2,  5,  10,  and  13  there  are  jumps  in  the  observation 
accuracy  for  heights.  The  jump  in  expected  error  at  abscissa  9 
is  due  to  the  complete  loss  wind  observations,  and  if  this 
abscissa  is  omitted,  then  the  jump  would  be  due  to  the  height 
accuracy  change  at  abscissa  10.  From  this  graph  it  is  apparent 
that  the  accuracy  of  height  observations  are  of  premier 
importance  in  the  height  analysis,  while  the  accuracy  of  the  wind 
observations  are  less  important  in  the  nominal  range  considered 
here,  except  that  not  having  wind  observations  at  all  also 
results  in  a  significant  increase  in  01  error. 

Significant  increases  in  error  of  the  analyzed  height  field 
are  seen  to  occur  with  decreasing  correlation  distance 
(increasing  value  of  a)  and  decreasing  constant  A,  for  all  three 
grid  and  observation  point  sets,  of  course,  the  effects  of 
observation  accuracy  are  considerably  smaller  for  the  sparse  NA 
grid.  For  all  three  grids  it  is  apparent  that  01  errors  are  more 
sensitive  to  changes  in  the  correlation  distance  (parameter  a) 
than  in  the  additive  constant  A. 
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abscissa 

h  error 

w  error 

1 

0 

1.0 

2 

5 

0.5 

3 

5 

1.0 

4 

5 

2.0 

5 

10 

0.0 

6 

10 

0.5 

7 

10 

1.0 

8 

10 

2.0 

9 

10 

00 

10 

20 

0.5 

11 

20 

1.0 

12 

20 

2.0 

13 

00 

1.0 

Table  3:  Abscissa  key  for  Figures  5-16 

The  corresponding  plots  for  the  errors  in  the  analyzed  wind 
fields  are  shown  in  Figures  8-10.  The  errors  shown  are  the  rms 
of  the  two  winds,  the  values  generally  being  quite  close 
together,  as  can  be  seen  by  referring  to  Tables  l.g.m  in  the 
Appendix  2.  The  errors  in  the  analyzed  wind  fields  generally 
follow  the  same  trend  as  the  errors  in  the  analyzed  height 
fields,  except  that  the  wind  errors  are  more  sensitive  to  the 
wind  observation  errors,  this  being  especially  prominent  when 
there  are  no  wind  observations  (at  abscissa  9).  One  interesting 
thing  to  notice  is  that  the  curves  for  correlation  functions  2-4 
are  very  close  together,  indicating  an  even  greater  dependence  on 
the  decay  rate  parameter,  a,  than  the  analyzed  height  errors 
exhibited. 

In  order  to  assess  the  importance  of  a  single  observation 
and  the  effects  of  many  missing  observations,  some  01  analysis 
errors  were  computed  based  on  missing  observed  values.  The 
missing  observations  are  shown  in  Figures  1-3  for  the  three 
grids.  The  expected  01  height  analysis  errors  are  shown  in 
Figures  11-13  for  the  three  grids,  and  the  expected  01  wind 
analysis  errors  in  Figures  14-16.  I  will  discuss  Figure  11  in 
some  detail,  and  a  similar  analysis  follows  for  the  other 
figures. 

The  figure  is  slightly  busy.  The  open  circles  and  open 
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squares  denote  the  values  for  correlation  functions  #4  and  #1, 
respectively,  the  same  values  as  given  in  Figure  5.  The  same 
symbol  with  the  +  denotes  the  analyzed  values  under  the  same 
assumptions  but  with  observed  values  missing  (appearing  twice; 
once  with  observation  of  winds  missing  at  the  location  noted 
previously,  the  second  with  observations  of  both  winds  and  height 
missing).  It  can  be  noted  that  the  missing  wind  observation  is 
nearly  undetectable  while  there  is  only  a  slight  degradation  with 
the  height  observation  also  missing.  This  is,  however,  more 
prominent  in  the  case  of  smaller  correlation  in  the  upper  graph 
for  correlation  function  #1.  The  shaded  symbols  denote  the 
values  obtained  when  one-half  (18)  of  the  observations  are 
missing  (again,  twice,  once  for  missing  winds,  and  once  for 
missing  heights  and  winds).  At  abscissa  9  the  symbol  overlays 
the  others  for  which  there  are  no  wind  observations  and  height 
observations  at  all  locations.  The  nominal  case  (abscissa  7) 
shows  about  a  15  percent  increase  in  the  analysis  errors  when  18 
wind  observations  are  missing,  and  about  a  50  percent  increase  in 
the  analysis  errors  for  each  correlation  function  when  the  entire 
observation  is  missing  at  18  points.  In  this  case,  a  15  percent 
increase  in  the  analysis  errors  corresponds  to  about  a  6  percent 
decrease  in  "skill”,  where  "skill"  is  taken  as  (l  -  expected 
analysis  error  relative  to  background  error),  and  the  50  percent 
increase  in  error  to  about  a  20  percent  decrease  in  skill.  The 
general  pattern  of  error  for  the  various  parameter  values  is 
generally  the  same  as  for  that  obtained  for  the  entire 
observation  set,  the  primary  difference  being  for  abscissa  9 
where  the  wind  observations  are  all  missing. 

The  Ol  analyzed  height  error  for  the  EC  grid  (Figure  12) 
follows  much  the  same  pattern,  except  that  the  increase  due  to 
one  missing  observation  is  somewhat  more  significant,  the  total 
number  of  observations  being  25  instead  of  36  as  in  the  MUS  grid. 
Again,  with  one-half  the  observations  missing,  the  error  is 
increased  by  about  50  percent,  and  the  drop  in  skill  about  25 
percent . 

The  MA  grid  results  in  Figure  13  differ  since  one  missing 
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observation  is  nearly  one-half  of  the  total  of  three.  Thus  there 
are  no  shaded  symbols  in  that  figure,  and  a  completely  missing 
observation  results  in  an  increase  of  height  analysis  error  of  a 
few  percent,  but  this  again  being  about  a  20  percent  drop  in 
skill . 

The  plots  for  the  wind  analysis  errors  are  given  in  Figures 
14-16.  Note  that  some  of  symbols  for  separate  correlation 
functions  are  overlaid  in  Figures  14  and  15.  Of  course,  the  wind 
analysis  is  sensitive  to  the  loss  of  one  or  more  observations, 
with  the  increase  in  wind  analysis  errors  in  Figure  14  for  the 
MUS  grid  showing  an  increase  of  about  6-8  percent  with  18  missing 
wind  observat j.ons ,  and  about  12-15  percent  when  the  entire 
observation  is  missing  at  18  points.  The  decrease  in  skill  here 
corresponds  to  about  20  percent  and  40  percent,  respectively.  On 
a  relative  basis  the  analysis  errors  are  significantly  larger  for 
winds,  with  smaller  relative  increases  in  the  error  when 
observations  are  missing,  however,  in  terms  of  skill  level,  the 
winds  are  more  dependent  on  the  observations.  Again,  the  general 
character  of  the  errors  follows  the  same  pattern,  with  the 
exception  of  abscissa  9,  corresponding  to  no  wind  observations. 

The  plots  for  the  wind  analysis  errors  for  the  EC  grid  in 
Figure  15  and  the  MA  grid  in  Figure  16  reveal  no  surprises.  The 
general  pattern  of  Figure  15  is  similar  to  Figure  14,  while 
Figure  16  reveals  that  a  very  small  skill  is  involved  in  this 
case,  so  missing  observations  have  little  affect. 

Unfortunately,  none  of  the  parameters  varied  above  are 
really  at  the  disposal  of  the  practitioner.  Still,  the  above 
Information  is  a  useful  aid  to  understanding  the  01  (and  SI) 
process  and  how  achievable  accuracy  is  affected  by  the  parameters 
in  the  process. 

The  more  important  practical  information  is  that  given  in 
Figures  17-22.  Once  again,  the  analysis  errors  are  plotted 
versus  a  single  abscissa  which  corresponds  to  various 
combinations  of  assumed  parameter  values.  The  nominal  background 
rms  error  is  30  m  for  the  height,  as  noted  in  the  previous 
section,  with  the  background  wind  errors  depending  on  the 
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correlation  function.  Nominal  height  observation  accuracy  is  10 
m,  and  nominal  wind  observation  accuracy  is  1  m/sec.  For  a  given 
spatial  correlation  function  number,  the  assumed  obseirvation 
accuracies  and  spatial  correlation  function  are  varied. 
Observation  accuracies  vary  by  a  factor  of  .5  or  2  from  the 
nominal,  while  correlation  function  number  varies  by  at  most  one. 

Consider  Figure  17.  The  three  "curves"  for  the  MUS  grid  and 
spatial  correlation  functions  2,  3,  and  4  will  be  discussed.  As 
in  previous  figures,  only  the  integer  abscissa  values  have 
meaning,  and  the  points  are  connected  only  to  enable  one  to  more 
easily  see  the  effects  of  changing  parameter  values.  It  is 
immediately  apparent  that  the  parameters  to  which  the  SI  scheme 
are  most  sensitive  are  embodied  in  abscissae  9,  12,  13,  and  to  a 
lesser  extent,  2,  4,  and  5.  Table  4  shows  the  relationship 
between  the  abscissae  and  the  parameter  variations,  and  we  see 
that  each  of  these  abscissae  except  9  are  for  assumed  height 
observation  errors  that  are  twice  the  nominal  value,  and 
abscissae  9,  12,  and  13  are  for  a  misspecified  correlation 
function  (greater  correlation).  Abscissa  -8,  -9,  -10,  -11,  -12, 
and  -13  also  show  relatively  larger  SI  errors,  and  each  of  these 
abscissa  are  for  low  assumed  correlation  as  well  as  improper 
assumed  height  observation  error.  It  appears  that  it  is  better 
to  underestimate  height  observation  errors  than  to  overestimate 
them,  although  there  is  a  peak  (but  smaller)  at  abscissa  -2, 
where  the  height  observation  error  is  underestimated. 

Note  that  the  graphs  for  correlation  functions  3  and  4  are 
quite  similar,  while  that  for  correlation  function  2  differs 
somewhat  for  negative  abscissae.  Table  1  shows  that  the  assumed 
correlation  function  (that  is  function  1)  for  large  negative 
abscissae  for  correlation  function  2  has  a  different  decay  rate, 
while  for  correlation  functions  3  and  4,  the  decay  rate  is  the 
same  as  that  of  2  and  3.  The  relatively  larger  effect  of  the 
improper  decay  rate  for  the  assumed  correlation  function  is  also 
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Table  4:  Abscissa  key  for  Figures  17-22 

apparent  in  the  SI  error  for  correlation  function  4  at  abscissae 
9,  12,  and  13,  as  noted  above.  Thus,  the  decay  rate  for  the 
correlation  function  seems  to  be  more  important  than  the  additive 
constant . 

For  the  EC  grid,  the  results  shown  in  Figure  18  indicate 
that  the  character  of  the  three  graphs  is  much  the  same.  The 
outstanding  difference  is  the  significantly  larger  SI  errors 
occurring  for  correlation  function  #2  when  the  assumed 
correlation  function  is  #1  (abscissae  -3  and  -6  to  -13).  Again, 
the  graph  for  correlation  function  #4  shows  larger  SI  errors  for 
assumed  correlation  function  #5  (larger  abscissae).  Both  of 
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these  cases  correspond  to  misspecif ied  decay  rates.  Looking 
higher  to  the  MA  grid,  we  see  same  effects:  cases  where  the 
decay  rate  for  the  correlation  function  is  misspecif ied  yield 
larger  increases  in  the  SI  height  error  than  when  the  additive 
constant  is  misspecif ied. 

The  results  for  the  MA  grid  shown  in  Figure  19  imply  that 
the  most  crucial  parameter  to  have  correct  in  such  sparse  regions 
is  the  decay  rate  for  the  correlation  function,  unfortunately  the 
most  difficult  to  estimate  in  such  cases.  Further,  it  appears  it 
is  probably  best  to  underestimate  the  decay  rate  in  sparse  (or 
semi-sparse  regions,  such  as  EC)  regions.  In  data  dense  regions 
the  height  errors  are  generally  less  sensitive  to 
misspecif ication  of  the  correlation  function  (with  the  exception 
of  this  being  in  combination  with  overestimates  of  the  height 
observation  error). 

The  corresponding  plots  for  SI  wind  errors  are  in  Figures 
20-22.  Here  the  behavior  of  the  errors  seems  to  be  less 
structured,  with  the  smallest  error  often  occurring  for  the 
correlation  function  corresponding  to  the  least  spatial 
correlation  (SI  height  errors  are  generally  a  decreasing  function 
of  correlation  function  #).  The  behavior  of  the  SI  wind  errors 
are  also  sensitive  to  the  misspecification  of  the  decay  rate  for 
the  correlation  function,  as  can  be  noted  by  the  correlation 
function  #2  values  for  large  negative  abscissae  and  for 
correlation  function  #4  for  large  positive  abscissae.  With  that 
exception,  the  general  behavior  of  the  SI  wind  errors  is  much  the 
same  as  for  the  SI  height  errors,  with  the  primary  dependence 
again  being  on  the  correct  specification  of  the  assumed  height 
error . 

One  additional  bit  of  information  can  be  squeezed  from  the 
data  generated  by  this  study.  This  concerns  the  relationship 
between  the  expected  error  based  on  the  01  calculations  versus 
the  actual  expected  error.  Of  course,  since  parameters  are 
estimated,  only  the  expected  error  for  01  can  be  calculated  when 
expected  errors  are  needed  (e.g.,  see  Goerss  (1989)).  We  give 
three  examples  to  show  how  this  proceeds,  and  to  show  the 


variation  in  the  values.  While  the  information  can  also  be 
obtained  by  looking  at  the  appropriate  figures,  the  information 
is  more  precisely  and  as  easily  obtained  from  the  tables  in 
Appendix  1. 

(1)  Let  us  suppose  that  the  assumed  values  of  the 
observation  errors  are  the  twice  the  nominal  ones  and  that  the 
assumed  spatial  correlation  function  is  #4.  Height  and  wind 
observation  errors  equal  20  m  and  2  m/sec,  respectively.  For 
the  MUS  grid  the  expected  error  for  the  heights  (from  Table 

1.1.4,  or  Figure  5,  abscissa  12  on  Cor  Ftn  4)  is  0.3462.  If 
the  actual  values  for  the  observation  errors  are  10  m  for 
heights  and  1  m  for  winds,  with  spatial  correlation  function 
#3,  then  the  actual  expected  error  for  the  MUS  grid  is  (from 
Table  6. 1.3-4,  or  Figure  17,  abscissa  13  on  Cor  Ftn  3)  is 
0.3092,  significantly  smaller  than  the  01  calculation  would 
indicate . 

(2)  Now  suppose  the  assumed  values  of  the  observation 
errors  are  one-half  the  nominal  values  and  that  the  assumed 
correlation  function  is  #3.  Height  and  wind  observation  errors 
equal  5  m  and  0.5  m/sec,  respectively.  For  the  MUS  grid  the 
expected  error  for  the  heights  (from  Table  1.1.3,  or  Figure  5, 
abscissa  2  on  Cor  Ftn  3)  is  0.2291.  Again  if  the  correct 
values  are  the  nominal  values  of  10  m  and  1  m/sec  for  the 
height  and  wind  observation  errors,  and  the  actual  spatial 
correlation  function  is  #4,  then  (from  Table  6. 1.4-3,  or  Figure 
17,  abscissa  -13  on  Cor  Ftn  4  )  the  correct  expected  error  is 
0.2819.  In  this  case  the  expected  errors  are  significantly 
larger  than  the  01  calculation  indicates. 

(3)  In  this  case,  suppose  the  spatial  correlation  function 
is  correct,  #4,  and  the  assumed  values  of  the  observation 
errors  are  5  m  and  2  m/sec  for  heights  and  winds,  respectively. 
Then  the  for  the  EC  grid  the  expected  height  error  (Table 

1.2.4,  or  Figure  6,  abscissa  4  on  Cor  Ftn  4)  is  0.3486.  If  the 
actual  observation  errors  are  10  m  and  l  m/sec  for  heights  and 
winds,  respectively,  then  the  actual  expected  error  (Table 

6. 2. 4-4,  or  Figure  18,  abscissa  -4  on  Cor  Ftn  4)  is  0.3870; 
again  the  01  expected  error  is  smaller. 

As  a  general  rule,  expected  error  as  calculated  by  01  will 
be  optimistic  when  the  observation  errors  are  underestimated  or 
when  the  spatial  correlation  is  overestimated.  The  latter 
usually  has  a  greater  Influence.  On  the  other  hand,  when 
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observation  errors  are  overestimated  or  when  the  spatial 
correlation  is  underestimated,  the  expected  error  estimates 
computed  by  01  will  be  pessimistic.  Any  expected  error  estimates 
from  operational  01  sources  should  be  treated  with  some  caution, 
the  examples  above  merely  serving  as  an  indication  of  the 
difficulties  and  not  as  a  guide  to  the  magnitude  of  the 
difference  between  the  computed  and  actual  values  that  may  occur 
in  practice. 

5.0  CONCLUDING  REMARKS 

The  results  of  this  study  demonstrate  that  SI  analyzed 
height  errors  are  more  sensitive  to  the  decay  rate  for  the 
spatial  correlation  function  than  for  the  additive  constant. 

While  this  study  concentrated  on  the  SOAR  correlation  function, 
similar  results  can  be  expected  for  other  correlation  functions 
which  are  controlled  by  parameters  governing  similar  properties. 
The  wind  errors  are  even  more  sensitive  to  proper  values  for  the 
decay  rate,  unsurprising  since  the  wind  correction  is  related  to 
the  derivative  of  the  spatial  correlation  function  for  the 
heights . 

Another  interesting  observation  is  that  it  is  better  to 
"make  the  same  mistake"  relative  to  the  observation  error  for  the 
heights  and  winds.  For  example,  if  observed  height  error  is 
specified  as  too  large,  smaller  analysis  errors  occur  when  the 
observed  wind  errors  are  also  specified  as  too  large,  rnther  than 
correct  or  too  small.  As  a  general  rule,  erring  on  the  side  of 
underestimating  the  observation  error  seems  to  result  in  smaller 
analysis  errors  than  erring  on  the  side  of  overestimation  of  the 
observation  error. 

The  effect  of  one  missing  wind  observation  is  vanishingly 
small.  This  is  true  even  when  one  observation  constitutes  a 
significant  portion  of  the  total  amount  of  data,  as  in  the  MA 
grid  case.  However,  two  things  come  into  play  in  this  case  to 
make  the  missing  data  still  rather  insignificant:  (1)  the 
missing  data  is  close  to  another  observation,  and  (2)  the  skill 
in  this  case  is  rather  low  anyway.  Missing  much  data  (about 
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one-half)  shows  significant  decreases  in  skill,  about  20-25 
percent  for  heights  and  up  to  40  percent  for  the  winds. 

The  relationship  between  the  expected  error  that  can  be 
computed  using  SI  and  the  actual  values  of  the  expected  error 
were  explored  briefly.  Practitioners  need  to  be  cognizant  of  the 
fact  that  these  two  values  may  be  significantly  different  from 
each  other. 


REFERENCES 

R.  Balgovind,  A.  Dalcher,  M.  Ghil,  and  E.  Kalnay  (1983),  A 
stochastic-dynamic  model  for  the  spectral  structure  of  forecast 
error  statistics,  Ifon.  Wea.  Rev.  Ill,  701-722. 

K.  H.  Bergman  (1979),  Multivariate  analysis  of  temperatures  and 
winds  using  optimum  interpolation.  Non.  Wea.  Rev.  107,1423-1444. 

R.  Franke  (1985),  Sources  of  error  in  objective  analysis.  Non. 
Wea.  Rev.  113,  260-270. 

R.  Franke  (1988),  Statistical  interpolation  by  iteration,  Mon. 
Wea.  Rev.  116,961-963. 

R.  Franke,  E.  Barker,  and  J.  Goerss  (1988),  The  use  of  observed 
data  for  the  initial-value  problem  in  numerical  weather 
prediction,  Coaput.  Nath.  Applic.  16,169-184. 

L.  S.  Gandin  (1965),  Objective  analysis  of  aeteorological  fields 
(1963),  translated  from  Russian  by  the  Israel  Program  for 
Scientific  Translations,  NTIS  TT65-50007. 

J.  Goerss  (1989),  Analysis  uncertainties  determination:  Navy 
Operational  Global  Atmospheric  Prediction  System  (NOGAPS) ,  TR# 
89-09,  Naval  Environmental  Prediction  Research  Facility, 
Monterey,  CA  93943. 

P.  Lonnberg  (1982),  Structure  functions  and  their  implications 
for  higher  resolution  analysis,  in  Proceedings  of  a  Workshop  on 
Current  Probleas  in  Data  Assimilation,  ECMHF,  pp.  142-178. 

A.  C.  Lorenc  (1981),  A  global  three  dimensional  multivariate 
st&tistical  interpolation  scheme,  Mon.  Wea.  Rev.  109,  701-721. 

A.  C.  Lorenc  and  0.  Mammon  (1988),  Objective  quality  control  of 
observations  using  Bayesian  methods.  Theory,  and  a  practical 
implementation,  Q.  J.  R.  Neteorol.  Soc.  114,  515-543. 

T.  W.  Schlatter  (1975),  Some  experiments  with  a  multivariate 
statistical  objective  analysis  scheme.  Non.  Wea.  Rev.  103,  246- 
257. 


17 


T.  W.  Schlatter,  G.  W.  Branstator,  and  L.  G.  Thiel  (1976), 

Testing  a  global  multivariate  statistical  objective  analysis 
scheme  with  observed  data.  Non.  Wea.  Rev.  104,  765-783. 

R.  S.  Seaman  (1983),  Objective  analysis  accuracies  of  statistical 
Interpolation  and  successive  correction  schemes,  Aust.  Net.  Nag. 
31,  225-240. 

H.  J.  Thi6baux,  (1985),  On  approximations  to  geopotential  and 
wind-field  correlation  structures,  Tellus  37A,  126-131. 

H.  J.  Thi6baux,  H.  L.  Mitchell,  and  D.  W.  Shantz  (1986), 
Horizontal  structure  of  hemispheric  forecast  error  correlations 
for  geopotential  and  temperature.  Non.  Wea.  Rev.  114,  1048-1066. 

H.  J.  Thi6baux  and  M.  A.  Redder  (1987),  Spatial  Objective 
Analysis,  Academic  Press,  London. 


18 


Figure  1:  The  Middle  US  (MUS)  grid  and  observation  locations. 

The  +'s  show  the  2.5°  grid.  The  open  squares  and  circles  and  the 
shaded  circle  show  observation  locations.  The  shaded  circle 
is  a  "missing  observation"  in  one  test  run,  and  all  circles  are 
"missing  observations"  in  another. 


Figure  2:  The  East  Coast  (EC)  grid  and  observation  locations. 
Symbols  as  in  Figure  1. 
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Figure  3:  The  Middle  Atlantic  (MA)  grid  and  observation 
locations.  Symbols  as  in  Figure  1,  except  there  is  only  a  single 
"missing  observation" . 


SPRTIRL  CORRELRTION 


Figure  4;  Spatial  correlation  functions. 


01  HEIGHT  ERRORS,  MR  GRID 


4  5  6  7  8  9  10  11  12  I 

PflRRMETER  VRRIRTIONS 

LEGEND 
o  -  COR  FT 
□  -  COR  FT 
^  -,COR  FT 
•  -  COR  FT 
■  -  COR  FT 


Figure  7:  As  in  Figure  5,  for  the  HA  grid  and 
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Figure  10:  As  in  Figure  8,  for  the  MA  grid  and  obser\ 
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01  H  ERROR  W/  MISSING  OBS,  MUS 
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Figure  11:  Expected  error  in  Ol  height  analysis  under  various 
observation  error  and  spatial  correlation  function  conditions  for 
the  MUS  grid  and  observation  set.  The  open  square  and  circle 
nark  the  values  for  the  entire  observation  set  for  correlation 
function  numbers  #1  and  #4,  respectively.  The  crossed  square  and 
circle  nark  two  values;  the  lower  for  one  missing  wind 
observation,  the  upper  one  for  both  height  and  wind.  The  shaded 
square  and  circle  mark  two  values;  the  lower  one  for  18  missing 
wind  observations,  the  upper  one  for  both  heights  and  winds. 
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Figure  13:  As  in  Figure  11, 
missing  observation  (only). 
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and 


*1  eT:^^ 


01  N  ERROR  W/  MISSING  OBS,  MUS 


PARAMETER  VARIATIONS 

LEGEND 

D  -  COR  FTN  1 
B  -  ONE  MISS 
■  -  18  MISS 
o  -  COR  FTN  4 
®  -  pNE  MISS 
•  -  18  MISS 

Figure  14:  Expected  error  in  Oi  wind  analysis  under  various 
observation  error  and  spatial  correlation  function  conditions  for 
the  MUS  grid  and  observation  set.  Symbols  as  in  Figure  11. 
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HEIGHT  SENSITIVITY  PLOTS,  MUS  GRID 


HEIGHT  SENSITIVITY  PLOTS,  EC  GRID 
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Figure  18:  As  in  Figure  17  for  the  EC  grid  and  observation 


HEIGHT  SENSITIVITY  PLOTS,  MR  GRID 
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Figure  19:  As  In  Figure  17  for  the  HA  grid  and  observation 
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20:  Expected  error  In  the  SI  wind  analysis  under  various 
correlation  function  and  parameter  values  (see  Table  4 
cissae  meanings),  HUS  grid  and  observation  set. 


NIND  SENSITIVITY  PLOTS,  EC  GRID 
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Figure  21:  As  in  Figure  20  for  the  EC  grid  and  observation 


WIND  SENSITIVITY  PLOTS,  Mfl  GRID 


Figure  22:  As  in  Figure  20  for  the  HA  grid  and  observation 


APPENDIX  1 :  TABLES 

The  tables  giving  the  values  plotted  in  the  various  figures, 
and  some  additional  data  as  well,  are  given  here.  Most  of  the 
information  is  more  readily  accessible  in  the  figures,  however 
the  tables  are  given  for  completeness. 

The  numbering  scheme  for  the  tables  is  a  key  to  the  grid  and 
observation  set  and  the  correlation  function  number  used.  The 
experiment  whose  results  are  given  by  the  tables  follows: 

Table  l.g.m:  Expected  error  under  various  values  of  parameters 
in  01,  nominal  observation  data. 

g  -  the  grid  used,  1  for  MUS,  2  for  EC,  3  for  MA. 
m  -  the  spatial  correlation  function  number,  1-5. 

Table  2.g.m:  Expected  error  under  various  values  of  parameters 
in  01,  wind  observation  missing  at  one  point,  as  noted  in  text 
and  Figures  1-3.  g  and  m  as  in  Table  l.g.m,  except  m=l  or  4. 

Table  3.g.m:  Expected  error  under  various  values  of  parameters 
in  01,  height  and  wind  observation  missing  at  one  point,  as  noted 
in  text  and  Figures  1-3.  g  and  m  as  in  Table  l.g.m,  except  m=l 
or  4. 

Table  4.g.m:  Expected  error  under  various  values  of  parameters 
in  01,  wind  observations  missing  at  about  one-half  the 
observation  points,  as  noted  in  text  and  Figures  1-3.  g  and  m  as 
in  Table  l.g.m,  except  m=l  or  4. 

Table  S.g.m:  Expected  error  under  various  values  of  parameters 
in  01,  height  and  wind  observations  missing  at  about  one-half  the 
observation  points,  as  noted  in  text  and  Figures  1-3.  g  and  m  as 
in  Table  l.g.m,  except  m=l  or  4. 

Table  6.g.B-n:  Expected  error  under  various  values  of  parameters 
in  SI. 

g  -  the  grid  used,  1  for  MUS,  2  for  EC,  3  for  MA. 
m  -  the  true  spatial  correlation  function  number,  1-5. 
n  -  the  assumed  spatial  correlation  function  number,  1-5, 
but  not  differing  more  than  one  from  m. 
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Tabla  1.1.3:  ixpaotad  arror  for  obaarvation  aocuraolaa 
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VabXa  1.3. Si  Brpaotad  arror  Cor  obaarva^ion  aoouraolaa 
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•rror^  hsiobln 

winds 

0.0 

S.O 

10.0 

20.0 

IfoBsiob^ 

0.0 

0.2682 

0.7040 

0.7014 

0.8 

0.2167 

0.2686 

0.3414 

O. 6764 

0.7048 

0.7389 

0.6706 

0.7023 

0.7346 

1.0 

0.1464 

O . 2177 

0.2667 

0.3429 

0.7604 

0.6870 

0.6779 

0.7061 

0.7377 

0,7711 

0.6447 

0.6720 

0.7034 

0.7403 

0.7797 

2.0 

0.2212 

0.2704 

0.3444 

0.6434 

0.7120 

0.7444 

0.6772 

0.7094 

0.7467 

HoWlnds 

0.3427 

0.4230 

0.4098 

Toblo  3.1.4S  Kxpootod  orror  Bor  obsorvoBien  ocourociM 


Ob« .  •rror 
winds 

0.0 

obs. 

8.0 

srror,  bsiobbs 

10.0  20.0 

]M9i^b 

0.0 

0.8937 

0.4886 

0.4871 

o.s 

0.8702 

0.4444 
O . 4447 

0.8934 

0.4887 

0.4871 

OOO 

SS9 

xl'Jb 

99  4 

MH  <4 

.  1.0 

0.8606 

0.4399 

0.4481 

0.8704 

0.4480 

0.4449 

0.8941 

O. 4889 
0.4874 

0.6800 
0.4783 
0.4748  ' 

0.4372 

0.9134 

0.9096 

a.o 

ft  S^ 

•489 

«SS 

ooo 

0.8980 

0.4867 

0.4842 

0.6811 

0.4771 

0.4783 

Vbifinds 

0.4631 

0.9121 

0.9190 

TobXo  a. 2.1:  Bryotod  orrov  Bor  oboorvoBloti  aocuroeloo 


Obm .  •rror 
winds 


0.0 


MoBsl^bt 


obs.  snror^  bsl^hes 

S.O  10.0  20.0 


0.0 

0.3744 

0.7741 

0.7762 

o.s 

0.3474 

O . 7S74 
0.7S4S 

0.37W2 

0.7745 

0.7766 

0.4344 

0.4014 

0.7444 

1.0 

0.3304 

0.7414 

0.7467 

0.3442 

0 . 7ft4S 

0 . 7S96 

0.3402 

0.7747 

0.7774 

0.4397 

0.4033 

0.7994 

0.4047 

0.4297 

0.4264 

2.0 

0.3&13 

O . 7626 
0.7640 

0.3440 

0.7442 

0.7426 

0.4444 

0.4045 

O . 4052 

MoWinde 

1 

0.4544 

0.4643 

0.4644 

T4bls  2.2.4:  Bx|^«o4«d  «)rror  Bor  observation  accuraeios 


obe .  error 
winds 

0.0 

ObS. 

6.0 

error »  heists 

10.0  20.0 

MOBeiqht 

0.0 

0.4420 

0.4464 

0.4470 

o.s 

0.4342 

0.4454 

0.9463 

0.9421 

0.4465 

0.4470 

0 . 4525 
0.4441 
0.4444 

1.0 

0.4364 

0.9450 

0.9461 

0.4343 

0.9454 

0.4463 

HAO 

fte^ 

ees 

see 

666 

0.4525 

0.4441 

0.4444 

or*  A 

eoe 

ess 

see 

666 

2.0 

0.9344 
0.4454 
O . 4464 

0.4423 

0.4465 

0.4471 

0.4527 

0.4441 

0.9444 

ICOHlnd. 

0.9520 

0.9496 

0.4922 

Table  2.3.1:  Expected  error  Bor  observation  aocuraeiee 


Obs .  error 
winds 

0.0 

Obs. 

5.0 

error*  heights 

10.0  20.0 

voBeight 

0.0 

0.4143 

0.4752 

0.4772 

0.5 

O . 4014 
O . 4741 
O . 4765 

0.4144 

0.9753 

0.4773 

0.4515 

0.9743 

0.4745 

1.0 

0.7464 

0.4734 

0.4764 

0.4017 

0.4743 

0.4767 

ses 

OAK 

ess 

ooo 

0.4514 

0.4745 

0.4747 

0.9405 
.  0.9464 
0.4462 

2.0 

0.4024 

0.4744 

0.4774 

0.4160 

0.4761 

0.4742 

0.4530 

0.9742 

0.9404 

Motlinds 

0.4374 

0.4447 

0.4914 

Table  2.3.4s  Expected  error  Bor  obeervation  acouraoiee 
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Gbm.  orror 
wind* 

0.0 

Obo. 

S.O 

orror#  bol^bto 

10.0  20.0 

1 

0.0 

0.4570 

0.4045 

0.4043 

0.5 

•  H4 

0  40 
H  OO 

4 

odd 

0.4571 

0.4044 

0.4044 

0.5400 

0.4310 

0.4320 

1.0 

0.4042 

0.7444 

0.7443 

0.4202 

0.7044 

0.7004 

0.4574 

0.4040 

0.4047 

0.5404 

0.4312 

0.4323 

0.7402 

0.4772 

O . 4444 

2.0 

0.42X4 

0.7055 

0.7017 

0.4547 

0.4040 

0.4057 

0.5414 

0.4324 

0.4333 

Mowindo  1 

0.5534 

0.4044 

0.4442 

1 

_ 1 

T«bX«  3.1.1s  B«p*ct*d  «rror  Bor  oboorvobion  aoeuraclMi 


Obo.  orror 
tflndo 

0.0 

Obo. 

5.0 

orror#  bol^bba 

10.0  20.0 

MoBslqlit 

0.0 

0.2441 

0.7034 

0.7047 

0.5 

0.2207 

0.6771 

0.4744 

0.2445 

0.7043 

0.7052 

0.3434 

0.7344 

0.7411 

1.0 

0.1910 

0.4546 

0.4533 

0.2217 

0.4747 

0.4754 

0.2447 

0.7054 

0.7047 

0.3454 

0.7344 

0.7424 

0.7414 

0.7444 

0.7414 

2.0 

0.2253 

0.4444 

0.4411 

0.2734 

0.7120 

0.7123 

0.3504 

0.7435 

0.7440 

Mowindo 

0.3442 

0.4244 

0.4133 

TobXo  3.l.4t  Bvp«cBod  orror  Bor  obMrvatlon  ocouraciM 


ob«.  arror 
winds 

0.0 

ObS. 

5.0 

•rror,  bsl^bs 

10.0  20.0 

moBrnl^Z 

0.0 

0.4000 

0.4545 

0.4542 

0.5 

0.5771 

0.4444 

0.4512 

0.4001 

0.4544 

0.4543 

0.4544 

0.4744 

0.4754 

1.0 

0.5474 

0.4452 

0.4474 

0.5773 

O . 4444 

0.4514 

0.4003 

0.4544 

0.4545 

0.4547 

0.4744 

0.4754 

0.4372 

0.4134 

0.4044 

?-o 

0.5742 

0.4504 

0.4522 

0.4013 

0.4404 

0.4403 

0.4554 

0.4744 

0.4747 

1 

1 

0.4723 

0.4144 

0.4221 

Tobia  3.2.1s  BrpocBad  orror  Bor  oboorvablon  aoouraoloo 
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Ob« .  orror 
vlnda 

0.0 

Ob«. 

5.0 

orror,  boigbto 

10.0  20.0 

KoBol^bb 

0.0 

0.3416 

0.7405 

0.7774 

0.5 

0.3504 

0.7514 

0.7604 

0.3420 

0.7404 

0.7742 

0.4406 

0.4024 

0.7441 

1.0 

0.3344 

0.7474 

0.7444 

0.3514 
0.7625 
O . 7621 

0.3431 

0.7421 

0.7745 

0.4420 

0.4043 

0.4006 

0.4047 

0 . 4297 
0.4264 

3.0 

0.3550 

0.7667 

0.7665 

0.3470 

0.7467 

0.7443 

0.4472 

0.4046 

0.4060 

MoMlnda 

pop 

00  6 

4 

MO  H 
HO  H 

Tabl«  3.2.41  arror  Bor  oboorvatlon  aoouraelM 


oto».  orror 

Obo. 

orror,  bolobbo 

vindo 

0.0 

5.0 

10.0 

20.0 

MoBaight 

0.0 

OflO 

400 

400 

ooo 

doo 

0.5 

0.4543 

0.4564 

0.4643 

0.4447 

0.4402 

0.4414 

0.4492 

0.4446 

0.9404 

1.0 

0.4534 

0.4544 

0.4570 

0.4643 

0.4445 

0.4446 

0.9497 

0.4402 

0.9415 

0.4457 

0.4491 

0.4443 

0.4446 

O . 4404 

0.4445 

2.0 

0.9545 

0.9572 

0.4645 

0.4444 

0.4402 

0.4415 

0.9443 

0.4497 

0.4904 

•rowindo 

0.4644 

0.4445 

0.9951 

Tablo  3.3.11  Bxpocf^od  orror  for  oboorvation  aoouraoloB 


Obs .  arror 
winds 

0.0 

Ob4. 

5.0 

srror,  bsi^dvts 

10.0  20.0 

KoBwlghb 

0.0  1 

40^ 

•  44 

•  4^ 

^  44 

ddd 

0.5 

^4S 

^  44 

ddd 

0.7440 

0.4645 

0.4747 

0.4264 

0.4755 

0.4777 

1.0 

4  4f| 

4  4h 

^  44 

d  d  d 

0.7750 

0.4661 

0.4734 

0.7443 

0.4647 

0.4744 

O . 4267 

0 . 4757 
0.4774 

0.4405 

0.4464 

0.4462 

2.0 

0.7761 

0.4666 

0.9741 

0.7403 

0.4702 

0.4756 

0.4277 

0.4762 

0.4746 

Motflndd 

0.4044 

0.4745 

0.4440 

Tablo  3.3.41  Brpacf  d  orror  for  obMrvablon  aocuraoloe 
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Oba .  arror 
winda 

0.0 

oba. 

B.O 

arror/  balobta 

10.0  20.0 

MoBaiobt 

0.0 

0.4933 

0.4457 

0.4406 

o.s 

0.4540 

0.4321 

0.4253 

0.4933 

0.4457 

0 . 4407 

990 

a^^ 

aaa 

666 

1.0 

0.4377 

0.t2<3 

0.4147 

0.4542 

0.4322 

0.4255 

0.4935 

0.4459 

0.4404 

0.5440 

0.4735 

0.4721 

0.4774 

0.9305 

0.9399 

3.0 

OO  0 

aaa 

bwn 

awb 

oa  a 

0.4942 

0.4465 

0.4413 

0.5444 

0.4741 

0.4727 

Noulnda 

0.5434 

0.4911 

0.4404 

Tabl«  4.1.1:  Bxp*c^4d  arror  toie  obMrva^ion  aeeuraolM 


Oba .  arror 

^>a. 

arror.  baidbra 

vinda 

0.0 

5.0 

10.0 

20.0 

MoBaldbb 

0.0 

0.3014 

0.7570 

0.7534 

0.5 

0.2477 

0.3016 

0.3920 

0.7244 

0.7573 

0.7943 

0.7167 

0.7537 

0.4016 

1.0 

0.2167 

0.2442 

0.3022 

0.3929 

0.4412 

O . 7026 

0.7253 

0.7543 

0.7994 

0.4541 

0.6924 

0.7175 

0.7545 

0.4025 

0.4644 

2.0 

0.2500 

0.3045 

0.3962 

0.7245 

0.7617 

0.4035 

0.7202 

0.7574 

0.4060 

itowlnda 

0.3427 

0.4230 

0.4095 

Tabla  4.1*4:  Bxpactad  arvor  for  obaorvation  aoouraeiaa 


Oba .  arror 

Oba. 

arror/  bai^ta 

vinda 

0.0 

5.0 

10.0 

20.0 

boBalobb 

0.0 

0.6297 

0.4447 

0.4460 

0.5 

0.6034 

0.6297 

0.6915 

0.4731 

0.4447 

0.9076 

0.4764 

0.4460 

0.9053 

1.0 

0.5932 

0.6039 

0.6294 

0.6916 

0.9039 

0.4679 

0.4732 

0.4444 

0.9077 

0.9537 

0.4729 

0.4769 

0.4462 

0.9054 

0.9441 

^.o 

0.6043 

0.6303 

0.6922 

• 

0.4735 

0.4452 

0.9041 

O . 4774 

0.4466 

0.9059 

Mowinda 

0.6632 

0.9121 

0.9190 

Tabla  4.2.1:  Bicpac^od  orror  for  ebaomratlon  aoouraeloa 


Obm .  arror 
vlnda 

0.0 

Oba. 

5.0 

arror,  haiqhba 

10.0  20.0 

MoBai^bb 

0.0 

0.4211 

0.4237 

0.4191 

0.5 

0.3426 
0.7977 
0 . 7970 

0.4212 

0.4239 

0.4193 

0.4911 

0.4555 

0.4443 

1.0 

0.3617 

0.7796 

0.7426 

0.3430 

0.7942 

0.7976 

0.4217 

0.4245 

0.4200 

0.4914 

0.4562 

0.4491 

0.4703 

0.4974 

0.4497 

2.0 

0.3443 
0.4003 
O . 7999 

0.4234 

0.4264 

0.4226 

0.4945 

0.4590 

0.4522 

MoWlnda 

0.4544 

0.4643 

0.4644 

Tabl*  4.2.4:  B:(p*cb«d  Mror  Bor  obsorvAtlon  accuraei** 
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Tabl«  S.1.1:  Bxp«o«*d  «rror  for  ebsorva^ion  aecuraeios 


Obm .  orror 
win4a 

0.0 

Oba. 

5.0 

orror,  bol^bta 

10.0  20.0 

NoBalqbb 

0.0 

0.4067 

0.4070 

0.4035 

0.5 

^aa 

odd 

ooa 

rt 

ooo 

oo  o 

ddd 

0.4729 

0.4263 

0.4249 

1.0 

0.3595 

0.7434 

0.7749 

0.3756 
O . 7936 
0.7445 

0.4040 

0.4045 

0.4049 

0.4741 

0.4275 

0.4300 

0.4412 

0.4541 

0.4644 

2.0 

0.3749 
0.7975 
O . 7920 

0.4114 

0.4126 

0.4045 

«•# 

ddd 

MoWinda 

0.4447 

0.4949 

0.4462 

Tablo  5.1.4:  Bxpoct.od  orror  for  obobrvation  acauraoioo 


Oba .  arror 
vinda 

0,0 

^a. 

5.0 

arrar,  baigb^a 

10.0  20.0 

lloBal9b9 

0.0 

0.7397 

0.9175 

0.9164 

0.5 

0.7260 

0.9123 

0.9121 

0.7394 

0.9175 

0.9164 

0.7753 

0.9246 

0.9260 

1.0 

0.7207 

0.9103 

0.9103 

a 

^  aa 

d  d  d 

0.7400 

0.9177 

0.9166 

0.7755 

0.9247 

0.9261 

0.9039 

0.9537 

0.9441 

2.0 

0.7264 

0.9129 

0.9127 

0.7406 

0.9141 

0.9171 

0.7762 

0.9292 

0.9266 

' 

HoMlnda 

0.7946 

0.9523 

0.9554 

Tablo  5.2.1:  Brpac^o4  orror  for  obaorvaflon  acouraeioa 
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Ob« .  ttriror 
vlnda 


0.0 


lloBttlgbt 


Ote.  Mrer, 

S.O  10.0  20.0 


0.0 

0.5221 

0.9590 

0.9540 

0.5 

0.5022 

0 . 9465 

O . 9440 

0.5224 

0.9592 

0.9543 

0.5574 

0.9739 

0.9592 

1.0 

0.4931 

0.439« 

0.9342 

00  0 

::s 

0.5232 

0.9590 

0.9552 

0.5594 

0.9745 

0.9591 

0.9703 

0.9979 

0.9997 

2.0 

0.5054 

0.9499 

0.9477 

0.5252 

0.9519 

0.9593 

0.5722 

0.9779 

0.9727 

MoWinda 

pop 

9901 
MK  • 
H9  9 
aH9 

Tabl*  S.2.4J  Bxpaetad  arror  toie  obMrvation  aocuraoiM 


•rr«  winds 


ASSUBSd 

5.0 

^s.  srrer^ 
10.0 

bsights 

20.0 

0.4531 

0.4474 

0.4717 

O . 4052 

0.4022 

0.4114 

0.4015 

0.7942 

0.4047 

0.4531 

0.4474 

0.4716 

0.4052 

O . 4022 

0.4117 

0.4015 

0.7942 

0.4047 

0.4532 

0.4474 

0.4713 

0 . 4052 

0.4022 

0.4117 

0.4015 

0.7942 

0.4047 

mne,  viiMta 


issnssd 

5.0 

Obs.  srror, 
10.0 

bsid^Ms 

20.0 

0.3045 

0.2945 

0.3279 

0.7034 

0.6939 

O . 7127 

0.7030 

0.6923 

0.7139 

0.3047 

0.2944 

0.3273 

0.7039 

0.6939 

0.7126 

0.7031 

0.6922 

0.7137 

0.3059 

0.2944 

0.3255 

0.7045 

0.6941 

0.7123 

0.7037 

0.6924 

0.7134 

Tsbls  6.1.1— 1:  Variatilon  of  oacpoc^od  srror 


Tsblo  6.l.2‘>3r  vsrlstlon  of  oxpmatmo  srror 


Assnd  Ctom. 
OTV,  winds 


Assussd 

5.0 

Obs.  srror, 
10.0 

bsiQhts 

20.0 

0.4543 

0.4709 

0.5476 

0.4064 

0.4120 

0.4379 

0.4020 

0.4041 

0.4371 

0.4540 

0.4699 

0.5457 

0.4061 

0.4116 

0.4374 

0.4017 

0.4077 

0.4366 

0.4564 

0.4666 

0.5347 

0.4054 

0.4105 

0.4356 

0.4011 

0 . 4067 

0.4349 

isssd  Obs. 
srr,  winds 


fcSSUSSd 

5.0 

Obs.  srror, 
10.0 

tasights 

20.0 

0.2945 

0.2793 

0 . 3016 

0.7149 

0.6940 

0.7104 

0.7152 

0.6964 

0.7116 

O . 2947 

0.2793 

0.3012 

0.7150 

0.6940 

0.7107 

0.7153 

0.6964 

0.7115 

0.2954 

0.2796 

0 . 2999 

0.7156 

0.6942 

0.7104 

0.7159 

0.6971 

0.7112 

Tsbls  6.1.1— a:  Vsirlstlon  of  sxpsctsd  srror 


Tsbls  6.1.3-as  Vsrlsblon  of  sspsotsd  srror 


srr,  winds 


Assussd 

5.0 

Obs.  srror, 
10.0 

hSiqhtS 

20.0 

0.3200 

0.3096 

0.3242 

0.7226 

O . 7064 

0.7032 

0.7232 

0.7054 

0.7014 

0.3201 

O . 3097 

0.3243 

0.7227 

0.7069 

0.7033 

0.7233 

0.7055 

0.7015 

0.3205 

0.3101 

0.3247 

0.7230 

0.7073 

0.7036 

0.7235 

0.7054 

0.7014 

Asssd  Obs. 
srr,  win<te 


Assussd 

5.0 

Obs.  srror, 
10.0 

bsiQbts 

20.0 

0.2922 

0.2749 

0.3059 

0.7121 

0.6977 

0.7139 

0.7122 

0.6965 

0.7152 

0.2925 

0.2749 

0.3054 

0.7122 

0.6976 

0.7137 

0.7123 

0.6965 

0.7150 

0*2937 

0.2791 

0.3034 

0.7129 

0.6979 

0.7134 

0.7130 

0.6967 

0.7147 

Tsbls  6.1.2-ls  Vsristlon  of  srpsotsd  srror 


Tsbls  6. 1.3-3:  Vsristlon  of  srpsotsd  srror 


Assussd 

5.0 

Obs.  srrsr, 
10.0 

bsi^bts 

20.0 

0.3064 

0.2941 

0.3217 

0.7060 

0.6936 

0.7094 

0.7054 

0.6920 

0.7100 

0.3066 

..  0.2941 

0.3212 

0.7060 

/  0.6936 

0.7093 

0.7055 

0.6920 

0.7099 

0.3076 

0.2942 

0.3197 

0.7066 

0.6934 

0.7090 

0.7060 

0.6921 

0.7096 

Asssd  Obs. 
srr,  winds 


ussussd 

5.0 

Obs.  srror, 
10.0 

bsi^bts 

20.0 

0.2901 

0.2793 

0.3114 

O . 7096 

0.6979 

0.7170 

0.7095 

0.6964 

0.7144 

0.2904 

0.2792 

0.3112 

O . 7097 

0.6979 

0.7164 

0.7096 

0.6964 

0.7146 

0.2914 

0.2792 

0.3092 

0.7106 

0.6942 

0.7165 

0.7104 

0.6970 

0.7143 

Tsbls  d.l.a-a^  VsristioA  of  srpsotsd  srror 


Tsbls  6.1.3*41  vsristlon  of  srpsotsd  srror 


iaaiinarl 

5.0 

obo.  orror, 
XO.O 

boi^tato 

20.0 

0.2tX» 

0.2655 

0.2464 

0.72X3 

0.70X4 

0.7X5X 

O . 7224 

0.70X3 

0.7X66 

0.2422 

0.2656 

0.2464 

0.72X4 

0.70X4 

0.7X50 

0.7225 

0.70X2 

0.7X64 

0.2435 

0.2654 

0.2450 

0.7222 

O.  7021 

0-7147 

0.7232 

0.70X5 

0.7X6X 

laounaf 

5.0 

Obn.  •rror, 
10.0 

b9iTbt« 

20.0 

0.2X47 

0.2016 

0.2237 

0.6547 

0.63X6 

0.6467 

0.6640 

0.6339 

0.65X7 

0.2X94 

0.20X5 

0.2224 

0.6595 

0.63X5 

0.6461 

0.6647 

0.6334 

0.6SXX 

0.2242 

0.2024 

0.2X92 

0.6645 

0.633X 

0.6462 

0.6694 

0.6354 

0.65X1 

TabI*  6.X.4>3 


variation  of  axpoctod  orror 


Tablo  tf.x.b^St  Variation  of  oxpootod  orror 


Aoood  obo. 
orr«  vlndo 


Araiinai 

5.0 

Obn.  «rror, 
XO.O 

bal^tn 

20.0 

0.2796 

0.2653 

0.29X7 

0.7X42 

0.70X6 

0.7X4X 

0.7X90 

O . 70X0 

0.7200 

0.2799 

0.2653 

0.29X2 

0.7X44 

0.70X6 

0.7X79 

0.7X92 

0 . 70X0 

0.7X94 

0.24X4 

0.2655 

0.2494 

0,7X43 

0.70X9 

0.7176 

0.720X 

0.70X3 

0.7X95 

Aoood  Obo. 
orr«  %rlndo 


»aaiina6 

5.0 

Obo.  nrror, 
XO.O 

bai^tn 

20.0 

0.5939 

0.5907 

0.6057 

0.4574 

0.4546 

0.4626 

0.4564 

0.4546 

O . 4604 

0 . 5940 

0.5907 

0.6056 

0.4574 

0.4546 

0.4626 

0.4564 

0.4546 

0.4607 

0.594X 

0.5907 

0.6054 

0.4574 

0,4546 

0.4626 

0.4565 

0.4546 

0 . 4607 

Tablo  6.x. 4*4:  Variation  of  oxpoctoo  orror 


Tablo  6.2.X-‘Xs  Variation  of  orpoctoo  orror 


jUfwntf 

5.0 

Ob«.  «rror/ 
XO.O 

baiobtr 

20.0 

0.2693 

0.2777 

0 . 3320 

0.70SX 

0.7X05 

0.7369 

0.7046 

0.7X07 

0.7420 

0.2695 

O . 2760 

0.3247 

0.7050 

0.7097 

0.7354 

0.7046 

0.7X00 

0.7404 

0.27X6 

O . 2724 

0.3X90 

0 . 7072 

0.7092 

0.7341 

0.7067 

O . 7095 

0.7391 

Aoood  Obo. 
orr,  vlnOo 


Xainaad 

5.0 

obn.  arror, 

XO.O 

baightn 

20.0 

0.63X1 

0.6364 

0.64X7 

0.4631 

0.4641 

0.4444 

0.4627 

0.4663 

0.4432 

0.6304 

0.6356 

0.6797 

0.4624 

0.4677 

0.4479 

0.4625 

0.4659 

0.4427 

0.6279 

0.63X3 

0.6725 

0.4620 

0.4665 

0.4463 

0.46X7 

0.4647 

0.44X0 

Tablo  6.x.4'*Ss  variation  of  oapoctod  orror 


Tablo  6.2.X-2S  variation  of  oxpoctod  orror 


AaaunaJ 

5.0 

Oba.  arror, 
10.0 

bai^bta 

20.0 

0.2343 

0.2095 

0.2XX7 

0.7005 

0.643X 

0.6370 

0.7044 

0.6465 

0.6394 

0.2344 

0.2099 

0.21X7 

0.7009 

0.6434 

0.6371 

0.7044 

0.6464 

0.6399 

0.2404 

0.2X13 

0.2X20 

0.7029 

0.6444 

0.6377 

0.7X06 

0.644X 

0.6404 

Anaunad 

5.0 

Oba.  arror, 
XO.O 

bai^ta 

20.0 

0.4642 

0.4655 

0.4423 

0.4059 

0.79X3 

0.7465 

0.7957 

0.7457 

0.7426 

0.4643 

0.4656 

0.4424 

0.4060 

0.7913 

0.7465 

0 . 7957 

0.7454 

0.7426 

0.4649 

0.466X 

0.4430 

0.4063 

0.79X6 

0.7464 

0.796X 

0.7461 

0.7429 

Tablo  6.1.S-4S  Variation  of  oapoctod  orror 


Tablo  e.a.a-ls  variation  of  orpootod  orror 


Ammmd 
•rr^  vinds 


Assussd 

5.0 

Obs.  orror, 
10.0 

hoi^hts 

20.0 

0.4351 

0.4274 

0.4454 

0.7414 

0.7713 

0.7430 

0.7750 

0.7643 

0.7745 

0.4353 

0.4277 

0.4455 

0.7414 

0.7713 

0.7429 

0.7760 

0.7643 

0.7744 

0.4362 

0.4275 

0.4445 

0.7423 

0.7715 

0.7427 

0.7765 

0.7645 

0.7742 

AM»d  Ote. 
•rr,  vinds 


>a«u— d 

5.0 

Obs.  srror, 
10.0 

boidbts 

20.0 

0.4041 

0.4015 

0.4230 

0.7446 

0.7750 

0.7490 

0.7744 

0.7717 

0.7441 

0.4044 

0.4014 

0.4226 

0.7447 

0.7749 

0.7449 

0.7749 

0.7716 

0.7440 

0 . 4096 

0.4015 

0.4213 

0.7455 

0.7752 

0.7446 

0.7797 

0.7719 

0.7434 

Variatldn  of  oxpoe^od  orror 


AM«d  Obm. 
mrie,  vinds 


Assiissd 

5.0 

Obs.  srror, 
10.0 

bsiqbts 

20.0 

0.4344 

0.4247 

0.4499 

0.7400 

0.7717 

0.7456 

0.7744 

0.7647 

0.7409 

0.4350 

0.4247 

0.4495 

0.7401 

0.7717 

0.7455 

0.7749 

0.7647 

0.7404 

0.4360 

0.4247 

0.4443 

0.7407 

0.7719 

0.7453 

0.7754 

0.7649 

0.7406 

Aooad  Ob«. 
•rr,  vlndo 


Assunsd 

5.0 

^s.  arror, 
10.0 

haiqhts 

20.0 

0.3464 

0.3771 

0.3924 

0.7942 

0.7740 

0.7479 

0.7464 

0.7745 

0.7432 

0.3471 

0.3771 

0.3921 

0.7943 

0.7740 

0.7477 

0.7466 

0.7744 

0.7430 

0.3443 

0.3774 

0.3912 

0.7951 

0.7743 

0.7475 

0.7473 

0.7747 

0.7424 

Tablo  6.2.2->3  3  Variation  of  axpacbad  arror 


Tabla  6. 2. 4-3:  Variation  of  oxpac^od  orroc 


•nr,  vlnda 


ASSUJMTd 

5.0 

Obs.  mrror, 
10.0 

bsiqbts 

20.0 

0.4111 

0.4025 

0.4193 

0.7495 

0.7753 

0.7444 

0.7424 

0.7720 

0.7403 

0.4113 

0.4025 

0  *  >1 

0.7496 

0.7753 

C  47 

0.7425 

0.7720 

O.  02 

0.4123 

0.4024 

0.4143 

0.7902 

0.7755 

0.7445 

0.7430 

0.7722 

0.7400 

Assikl  Ob«. 
•rr,  winds 


ASSItfMKS 

5.0 

Obs.  srror, 
10.0 

bsiqbts 

20.0 

0.3454 

0.3767 

0.3944 

0.7914 

0.7774 

0.7499 

0.7444 

0.7742 

0.7451 

0.3457 

0.3767 

0.3944 

0.7916 

0.7774 

0.7494 

0.7446 

0.7742 

0.7449 

0.3470 

0.3769 

0.3933 

0.7925 

0.7740 

0.7496 

0.7454 

0.7745 

0.7447 

Tabl«  6. 2. 3-2:  Variation  of  •xp«ct«d  orror 


Tabl«  4. 2.4-4:  Variation  of  oxpootod  orror 


•rr ,  winds 


Assusad 

5.0 

Obs .  srror , 
10.0 

bsiqbts 

20.0 

0.4091 

0.4011 

0.4192 

0.7464 

0.7747 

0,7467 

0.7403 

0.7714 

0.7420 

0.4093 

0.4011 

0.4144 

0.7469 

/  0.7747 

0.7466 

0.7404 

0.7714 

0.7414 

0.4104 

0.4012 

0.4174 

0.7475 

0.7749 

0.7463 

0.7410 

0.7716 

0.7416 

ASSMd  obs. 

•rr,  winds 


Assuasd  Obs.  srror,  bsl^bts 


Tabl«  6. 2. 3-3:  Variation  of  •xp«et«d  orror 


TabXs  6.2.4-51  Variation  of  •xpsotsd  orror 


Ob«. 
•rr,  winds 


0.3104 

0.2946 

0.3036 

0.7703 

0.7240 

0.7203 

0.7549 

0.7196 

0.7167 

0.3112 

0.2950 

0.3034 

0.770S 

0.7243 

0.7204 

O.  7553 

0.7199 

0.7164 

0.3136 

0.2967 

0.3049 

0.7729 

0.7254 

0.7212 

0.7574 

0.7213 

0.7176 

Assuaad 

5.0 

Obs.  arror, 
10.0 

lialgtots 

20.0 

0.4449 

0.4913 

0.9013 

0.9674 

0.9669 

0.9670 

0.9740 

0.9732 

0.9723 

0.4490 

0.4914 

0.9014 

0.9674 

0.9669 

0.9671 

0.9740 

0.9732 

0.9724 

0.4492 

0.4916 

0.9016 

0.9679 

0.9670 

0.9672 

0.9741 

0.9733 

0.9724 

Tabls  6. 2. 5-4:  Variation  ot  axpacbad  arror 


Tania  6. 3. 2-1:  Variation  of  arpaetad  arror 


Aaaad  Obm. 
arr«  winds 


Assuaad 

5.0 

^s .  arror  , 
10.0 

haidhts 

20.0 

0.2952 

0.7322 

0.7253 

0.2444 

0.7119 

0.7047 

awn 
wn  a 
on  w 

odd 

0.2961 
0.7329 
O  .  7260 

0.2443 

0.7114 

0.7046 

0.3009 

0.7226 

0.7144 

0.3001 

0.7371 

0.7300 

0.2452 

0.7132 

0.7100 

O  .  2990 
0.7224 
0.7190 

arr,  winds 


Assuaad 

5.0 

M>s.  arror, 
10.0 

balqbts 

20.0 

0.4541 

0.4572 

O . 4620 

0.9595 

0.9547 

0.9603 

0.9653 

0.9651 

0.9654 

0.4541 

0.4572 

0.4619 

0.9595 

0.9547 

O  .  9603 

0.9653 

0.9651 

0.9654 

0.4543 

0.4573 

0.4620 

0.9596 

0.9547 

0.9603 

0.9653 

0.9651 

0.9654 

Tabla  6.2. 5-5:  variation  of  axpaetad  arror 


Tabla  5. 3.2-2:  Variation  of  axpaetad  arror 


Assad  Obs. 
arr,  winds 


Assuaad 

5.0 

Obs .  arror , 
10.0 

halqhts 

20.0 

0.9359 

0.9356 

0.9379 

0.9419 

0.9417 

0.9425 

0.9444 

0.9443 

0.9447 

0.9359 

0.9356 

0.9379 

0.9419 

0.9417 

0.9425 

0.9444 

0.9443 

0.9447 

0.9359 

0.9356 

0.9379 

0.9419 

0.9417 

0.9425 

0.9444 

0.9443 

0.9447 

Assad  Obs. 
arr^  winds 


iSBuaad 

5.0 

Obs.  arrar, 
10.0 

halqbts 

20.0 

0.4640 

0.4649 

0.4639 

0.9595 

0.9590 

0.9611 

0.9656 

0.9656 

0.9667 

0.4640 

0.4649 

0.4634 

0.9595 

0.9590 

0.96X1 

0.9656 

0.9656 

0.9667 

0.4641 

0.4649 

0.4639 

O  .  9596 

0.9590 

0.9611 

0.9657 

0.9656 

0.9664 

Tabla  6. 3. 1-1:  variation  of  axpaetad  arror 


Tabla  6. 3. 2-3:  Variation  of  axpaetad  arror 


Assad  OkM. 
arr,  winds 


Assuaad 

5.0 

Obs .  arror , 
10.0 

halqbts 

20.0 

0.9916 

0.9473 

0.9407 

0.9474 

0.9479 

0.9497 

0.9901 

0.9902 

0.9904 

0.9907 

0.9463 

0.9796 

0.9476 

0.9474 

0.9495 

0.9499 

•  0.9900 

0.9907 

0.9472 

0.9426 

0.9757 

0.9471 

0.9472 

0.9449 

0.9494 

0.9495 

0.9901 

Assad  Obs. 
arr,  winds 


Assxiaad 

5.0 

Obs.  arror, 
10.0 

balqhts 

20.0 

0.4244 

0.4294 

O . 4402 

0.9616 

0.9604 

0.9613 

0.9643 

0.9677 

0.9675 

0.4244 

0.4294 

O . 4402 

0.9616 

0.9604 

0.9613 

•  0.9643 

0.9677 

0.9675 

0.4246 

0.4295 

0.4402 

0.9617 

0.9604 

0.9613 

0.9643 

0.9677 

0.9676 

Tabla  s. 3.1-2:  variation  of  axpaetad  arror 


Tabla  6. 3. 3-2:  Variation  of  axpaetad  arror 


Aasad  Ob«. 
err,  vinds 


Aeeuaed 

5.0 

Obe.  error, 
10.0 

beigbte 

20.0 

0.4212 

0.4201 

0.4261 

0.9610 

0.9600 

0.9617 

0.9673 

0.9671 

0.9677 

0.4212 

0.4201 

0.4261 

0.9610 

0.9600 

0.9617 

0.9673 

0.9670 

O  .  9677 

0.4213 

0.4202 

0.4261 

0.9611 

0.9601 

0.9617 

0.9673 

0.9671 

0.9674 

Aeeuaed 

5.0 

Obe.  error, 
10.0 

bei^hte 

20.0 

0.4047 

0.4001 

0.7945 

0.9664 

0 . 9667 

0.9644 

0.9733 

0.9734 

0.9743 

0 . 4030 

0.7943 

0.7967 

0.9664 

0.9662 

0.9679 

0.9724 

0.9729 

0.9739 

0.7944 

0.7930 

0.7912 

0.9653 

0.9647 

0.9663 

0.9717 

0.9717 

0.9727 

Table  d.S.S'-d:  Variation  ot  expected  error 


Table  6.3.4*5t  Variation  of  expected  error 


Aeead  Obe. 
err,  vinde 


Aee^iaed 

5.0 

Obe.  error, 
10.0 

heidhte 

20.0 

0.4266 

0.4234 

0.4249 

0 . 9609 

0.9602 

0.9623 

0.9674 

0.9673 

0.9643 

0.4266 

0.4234 

0.4249 

0.9609 

0 . 9602 

0.9623 

O  .  9674 

0.9673 

0.9643 

0.4264 

0.4234 

0.4249 

0.9610 

O . 9602 

0.9623 

0.9675 

0.9674 

0.9644 

Aeead  obe. 
err,  winde 


Aeeuaed 

5.0 

^Mi.  error, 
10.0 

beidbte 

20.0 

0.6957 

0.6940 

0.7169 

0.9403 

0.9370 

0.9347 

0.9536 

0.9524 

0.9522 

O  .  6960 

0.6943 

0.7172 

0 . 9405 

0.9372 

0.9344 

0.9534 

0.9S2S 

0.9523 

0.6972 

0.6995 

'>.7144 

0.9412 

0.9370 

0.9394 

0.9544 

0.9531 

0.9529 

Table  6. 3. 3-4$  Variation  of  expected  error 


Table  6.3.5~4:  Variation  of  expected  error 


Aeead  obe. 
err,  wlnde 


Aeeuaed 

5.0 

Obe.  error, 
10.0 

heid^te 

20.0 

0.7410 

0.7414 

0.7940 

0.9626 

0.9612 

0.9624 

0.9692 

0.9647 

0.9649 

0.7410 

0.7414 

0.7940 

0,9626 

0.9612 

0.9624 

O . 9692 

0.9647 

0.9649 

0.7412 

0.7419 

0.7940 

0.9627 

0.9613 

0.9624 

O . 9692 

0.9647 

0.9649 

err,  winde 


ieeuaed 

5.0 

Obe.  error, 
10.0 

beidhte 

20.0 

0.6442 

0.6424 

0.6925 

0.9324 

0.9307 

0.9332 

0.9473 

0.9469 

0.9479 

0.6443 

0.6424 

0.6924 

0.9329 

0.9306 

0.9331 

0.9473 

0.9464 

0.9474 

0.6452 

0.6424 

0.6924 

0.9334 

0.9310 

0.9334 

0.9477 

0.9471 

0.9441 

Table  6.3. 4*3 


Variation  of  expected  error 


Table  6.3.5*58  variation  of  expected  error 


Bead  ^e. 

rr ,  vinde 


Aeeuaed  obe.  error,  lieidlite 


5.0 

10.0 

20.0 

0.7747 

0.7775 

0.7445 

0.9621 

0.9610 

0.9624 

0.9646 

0.9644 

0.9691 

0.7747 

0.7775 

0.7445 

0.9622 

O . 9610 

0.9627 

0.9646  . 

0.9643 

0.9691 

0.7749 

0.7776 

0.7446 

0.9623 

0.9611 

0.9624 

0.9647 

0.9644 

0.9691 

variation  of  expected  error 


oooooooon  oonooooooonnoonooonoooooo 


APPENDIX  2:  SUBROUTINE  TCOR 


This  appendix  contains  the  listing  for  the  subroutine  that 
evaluates  the  true  covariance  and  cross-covariance  function 
values  between  a  given  (grid)  point  and  the  observation  points. 


SUBROUTINE  TCOR ( NCF , THD , PHD , MH , TUD , PUD , NU , TVD , PVD , NV , TD , PD , KH , KU , 
1  KV,VAR,CAPPA,CVM,NCVM) 

THIS  SUBROUTINE  CALCULATES  THE  TRUE  COVARIANCE  BETWEEN  THE 
OBSERVATION  VARIABLES  AND  INDICATED  VARIABLES  AT  A  GIVEN  POINT 

THE  ARGUMENTS  ARE 


INPUT  ARGUMENTS 
NCF 

THD,PHD;NH 

TUD,PUD;NU 

TVD,PVD;NV 


KH 

KU 

KV 

VAR 

CAPPA 

TD,PD 

NCVM 


COVARIANCE  FUNCTION  NUMBER 

ARRAY  OF  HEIGHT  OBSERVATION  LOCATIONS,  NH  OF  THEM 

ARRAY  OF  U-WIND  OBSERVATION  LOCATIONS,  NU  OF  THEM 

ARRAY  OF  V-WIND  OBSERVATION  LOCATIONS,  NV  OF  THEM 

THE  ABOVE  ARE  IN  DEGREES,  LONGITUDE  AND  LATITUDE, 
RESPECTIVELY 

NONZERO  IF  H  COVARIANCE  TO  BE  COMPUTED 
NONZERO  IF  U  COVARIANCE  TO  BE  COMPUTED 
NONZERO  IF  V  COVARIANCE  TO  BE  COMPUTED 
VARIANCE  OF  HEIGHT-HEIGHT  ERRORS 
CORIOLIS  CONSTANT 
GIVEN  (GRID)  LOCATION,  DEGREES 
ROW  DIMENSION  OF  CVM  ARRAY 


OUTPUT  ARGUMENT 

CVM  -  ARRAY  OF  COVARIANCES  BETWEEN  OBS  VARS  AND  GRID  LOGS 

AN  NH+NU+NV  BY  3  ARRAY 

IMPLICIT  REAL* 8  (A-H,0-Z) 

PARAMETER  (NSZ=36) 

DIMENSION  THD(NH) ,PHD(NH) ,TUD(NU) ,PUD(NU) ,TVD(NV) ,PVD(NV) , 

1  CVM(NCVM,3) 

DIMENSION  TOH(NSZ) ,POH(NSZ) ,TOU(NSZ) ,POU(NSZ) ,TOV(NSZ) ,POV(NSZ) 
DIMENSION  AP(5),AAP(5) 

COMMON  /lO/KOUT 
DATA  NTIME/0/ 

DATA  AP , AAP/5D-6 , 3 . 0825D-6 , 3 . 0825D-6 , 3 . 0825D-6 , 2 . 188D-6 , 

1  ODO,ODO, .15D0, .2722D0, .2722D0/ 

THIS  ROUTINE  IS  SET  UP  TO  ACCEPT  GENERAL  (ISOTROPIC)  COVARIANCE 
FUNCTIONS.  THIS  REQUIRES  THE  DEFINITION  OF  A  NtJMBER  OF  ARITHMETIC 
STATEMENT  FUNCTIONS  WHICH  DEFINE  THE  DISTANCE  FUNCTION  IN  TERMS  OF 
LATITUDE  AND  LONGITUDE  (PHI  AND  THETA),  ITS  DERIVATIVES  WRT  PHI  AN 
THETA,  AS  WELL  AS  THE  COVARIANCE  FUNCTION  AND  IT'S  DERIVATIVES. 

THE  FUNCTIONS  NEEDED  ARE 
DIST  -  DISTANCE  FUNCTION 
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ooo  no  ooooooooooo 


cov 

DCOV 

D2C0V 

VARU 

VARV 

DSDP 

DSDT 

DDPl 

DDP2 

DDP3 


ISOTROPIC  COVARIANCE  FUNCTION  AS  A  FUNCTION  OF  DISTANCE 
DERIVATIVE  OF  COV  WRT  DISTANCE 

2ND  DERIVATIVE  OF  COV  WRT  DISTANCE  -  DCOV/DISTANCE 

VARIANCE  OF  U-WIND 

VARIANCE  OF  V-WIND 

PARTIAL  DIST  WRT  PHI I 

PARTIAL  DIST  WRT  THETAI 

PART  OF  2ND  PARTIAL  DIST  WRT  PHII  AND  PHIJ 
PART  OF  2ND  PARTIAL  DIST  WRT  THETAI  AND  THETAJ 
PART  OF  2ND  PARTIAL  DIST  WRT  PHII  AND  THETAJ 


DATA  RAD/6.372D6/ 

DIST ( P , PI , T , T1 ) =RAD*DSQRT ( ( P-Pl ) **2+ ( ( T-Tl ) *COS ( ( P+Pl ) /2 . DO ) ) **  2 ) 
COV(S)  =  VAR*(OMAA*(1.DO+A*S)*EXP(-A*S)+AA) 

DCOV(S)  =  -VAR*OMAA*A*A*S*EXP(-A*S) 

D2COV(S)  =  VAR*OMAA*A**3*S*EXP(-A*S) 

VARU(P)  =  VAR*OMAA*(RAD*A)**2 
VARV(P)  =  VAR*OMAA*(RAD*A*COS(P) )**2 

DSDP(P,P1,T,T1,S)  =  RAD**2/S*(-(T-T1)**2*SIN(P+P1)/4.D0+(P-P1) ) 
DSDT(P,P1,T,T1,S)  =  RAD**2/S*(T“Tl)*COS( (P+P1)/2.D0)**2 
DDP1(P,P1,T,T1,S)  =  RAD**2/S*(1.D0  +  (T-Tl ) **2*COS(P+Pl )/4 . DO ) 
DDP2(P,P1,T,T1,S)  =  (RAD*COS( (P+P1)/2.D0) )**2/S 
DDP3(P,P1,T,T1,S)  =  RAD**2/S*(T-T1)*SIN(P+P1)/2.D0 
IF(NTIME.NE,NCF)  THEN 

SET  PARAMETERS  FOR  THIS  FUNCTION  NUMBER 


NCFR=MAX(1,NCF) 

NCFR=MIN(NCFR,5) 

A  *  AP(NCFR) 

AA  =  AAP(NCFR) 

OMAA  =  l.DO-AA 
WRITE(KOUT,l)A,AA 
PI=DATAN ( 1 . DO ) *4 . DO 
DTR  =  PI/180. DO 
NTIME  =  NCF 

ENDIF 

CONVERT  LOCATIONS  TO  RADIANS 


DO  100  1=1, NH 

POH(I)  =  PHD(I)*DTR 
TOH(I)  =  THD(I)*DTR 
100  CONTINUE 

DO  110  1=1, NU 

POU(I)  =  PUD(I)*DTR 
TOU(I)  =  TUD(I)*DTR 
110  CONTINUE 

DO  120  1=1, NV 

POV(I)  =  PVD(I)*DTR 
TOV(I)  =  TVD(I)*DTR 
120  CONTINUE 
P  =  PD*DTR 
T  =  TD*DTR 
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C  START  CALCULATIONS 

K  =  0 

SP  =  SIN(P) 

CP  =  COS(P) 

IF(NH.NE.O)  THEN 
DO  200  1=1, NH 
K  =  K  +  1 

S  =  DIST(P,POH(I) ,T,TOH(I)) 

DF  =  DCOV(S) 

IF(KH.GT.O)  CVM(K,1)  =  COV(S) 

CVM(K,2)  =  O.DO 
CVM(K,3)  =  O.DO 
IF(KU.GT.O  .AND.  S.GT.O.DO) 

1  CVM(K,2)  =  DF*DSDP(P,POH(I) ,T,TOH(I) ,S) 

IF(KV.GT.O  .AND.  S.GT.O.DO) 

1  CVM(K,3)  =  DF*DSDT(P,POH(I) ,T,TOH(I) ,S) 

CVM(K,2)  =  -CVM(K,2)*CAPPA/SP 
CVM(K,3)  =  CVM(K,3)*CAPPA/CP/SP 
200  CONTINUE 
ENDIF 

IF(NU.NE.O)  THEN 
DO  220  1=1, NU 
K  =  K  +  1 

S  =  DIST(P,POU(I) ,T,TOU(I) ) 

DF  =  DCOV(S) 

DD2  =  D2COV(S) 

SPI  =  SIN(POU(I)) 

CPI  =  COS(POU(I)) 

CVM(K,1)  =  O.DO 
IF(KH.GT.O  .AND.  S.GT.O.DO) 

1  CVM(K,1)  =  DF*DSDP(POU(I),P,TOU(I),T,S) 

CVM(K,2)  =  VARU(P) 

IF(KU.GT.O  .AND.  S.GT.O.DO) 

1  CVM(K,2)  =  -DF*DDP1(P0U(I) ,P,TOU(I) ,T,S)  + 

2  DD2*DSDP{P,POU(I),T,TOU(I) ,S)*DSDP(POU(I) ,P,TOU(I) ,T,S) 
CVM(K,3)  =  O.DO 

IF(KV.GT.O  .AND.  S.GT.O.DO) 

1  CVM(K,3)  =  -DF*DDP3(P,POU(I) ,T,TOU(I) ,S)  + 

2  DD2*DSDP(POU(I) ,P,TOU(I) ,T,S)*DSDT(P,POU(I) ,T,TOU(I) ,S) 
CVM(K,1)  =  -CVM(K,1)*CAPPA/SPI 

CVM(K,2)  =  CVM(K,2)*CAPPA**2/SP/SPI 
CVM(K,3)  =  -CVM(K,3)*CAPPA**2/SP/CP/SPI 
220  CONTINUE 
ENDIF 

IF(NV.GT.O)  THEN 
DO  240  1=1, NV 
K  =  K  +  1 

S  =  DIST(P,POV(I),T,TOV(I)) 

DF  =  DCOV(S) 

DD2  =  D2COV(S) 

SPI  =  SIN(POV(I)) 

CPI  =  COS(POV(I)) 

CVM(K,1)  =  O.DO 
IF(KH.GT.O  .AND.  S.GT.O.DO) 
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CVM(K,1)  =  DP*DSDT(POV(I) ,P,TOV(I) ,T,S) 

CVM(K,2)  =  O.DO 
IF(KU.GT.O  .AND.  S.GT.O.DO) 

1  CVM(K,2)  =  -DF*DDP3(POV(I) ,P,TOV(I) ,T,S)  + 

2  DD2*DSDP(P,POV(I) ,T,TOV(I) ,S)*DSDT(POV(I) ,P,TOV(I) ,T,S) 
CVM(K,3)  =  VARV(P) 

IF(KV.GT.O  .AND.  S.GT.O.DO) 

1  CVM(K,3)  =  ~DF*DDP2(P,POV(I) ,T,TOV(I) ,S)  + 

2  DD2*DSDT(POV(I) ,P,TOV(I) ,T,S)*DSDT(P,POV(I) ,T,TOV(I) ,S) 
CVM(K,1)  =  CVM(K,1)*CAPPA/CPI/SPI 

CVM(K,2)  =  -CVM(K,2)*CAPPA**2/CPI/SPI/SP 
CVM(K,3)  =  CVM(K,3)*CAPPA**2/CPI/CP/SP/SPI 
240  CONTINUE 
ENDIF 
RETURN 

1  FORMAT('  TCOR  -  01/12/90:  2ND  AR,  A,AA  =  ',1P,2E11.3) 

END 
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INITIAL  DISTRIBUTION  LIST 


ASST.  FOR  ENV.  SCIENCES 
ASST.  SEC.  OF  THE  NAVY  (R&D) 
ROOM  5E731,  THE  PENTAGON 
WASHINGTON,  DC  20350 

OFFICE  OF  NAVAL  RESEARCH 
CODE  422AT 
ARLINGTON,  VA  22217 

COMMANDING  OFFICER 
FLENUMOCEANCEN 
MONTEREY,  CA  93943 

NAVAL  POSTGRADUATE  SCHOOL 
METEOROLOGY  DEPT. 

MONTEREY,  CA  93943 

LIBRARY  ( 2 ) 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 

COMMANDER 

NAVAIRSYSCOM  (AIR-330) 
WASHINGTON,  DC  20361 

USAFETAC/TS 

SCOTT  AFB,  IL  62225 

AFGWC/DAPL 

OFFUTT  AFB,  NE  68113 
AFGL/OPI 

HANSCOM  AFB,  MA  01731 

COMMANDER  &  DIRECTOR 
ATTN;  DELAS-D 
U.S.  ARMY  ATMOS.  SCI.  LAB 
WHITE  SAND  MISSILE  RANGE 
WHITE  SANDS,  NM  88002 

NOAA-NESDIS  LIAISON 
ATTN:  CODE  SC2 

NASA-JOHNSON  SPACE  CENTER 
HOUSTON,  TX  77058 

ACQUISITIONS  SECT.  IRDB-DB23 
LIBRARY  &  INFO.  SERV. ,  NOAA 
6009  EXECUTIVE  BLVD. 
ROCKVILLE,  MD  20852 


NATIONAL  WEATHER  SERVICE 
WORLD  WEATHER  BLDG.,  RM  307 
5200  AUTH  ROAD 
CAMP  SPRINGS,  MD  20023 

DIRECTOR 

GEOPHYS.  FLUID  DYNAMICS  LAB 
NOAA,  PRINCETON  UNIVERSITY 
P.O.  BOX  308 
PRINCETON,  NJ  08540 

LABORATORY  FOR  ATMOS.  SCI. 
NASA  GODDARD  SPACE  FLIGHT  CEN. 
GREENBELT,  MD  20771 

COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN;  DR.  WILLIAM  GRAY 
FORT  COLLINS,  CO  80523 

CHAIRMAN,  METEOROLOGY  DEPT. 
UNIVERSITY  OF  OKLAHOMA 
NORMAN,  OK  73069 

COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN;  LIBRARIAN 
FT.  COLLINS,  CO  80523 

UNIVERSITY  OF  WASHINGTON 
ATMOSPHERIC  SCIENCES  DEPT. 
SEATTLE,  WA  98195 

FLORIDA  STATE  UNIVERSITY 
ENVIRONMENTAL  SCIENCES  DEPT. 
TALLAHASSEE,  FL  32306 

DIRECTOR 

COASTAL  STUDIES  INSTITUTE 
LOUISIANA  STATE  UNIVERSITY 
ATTN;  0.  HUH 
BATON  ROUGE,  LA  70803 

UNIVERSITY  OF  MARYLAND 
METEOROLOGY  DEPT. 

COLLEGE  PARK,  MD  20742 

CHAIRMAN 

METEOROLOGY  DEPT. 
MASSACHUSETTS  INSTITUTE  OP 
TECHNOLOGY 
CAMBRIDGE,  MA  02139 
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ATMOSPHERIC  SCIENCES  CENTER 
DESERT  RESEARCH  INSTITUTE 
P.O.  BOX  60220 
RENO,  NV  89506 

CENTER  FOR  ENV.  &  MAN,  INC. 
RESEARCH  LIBRARY 
275  WINDSOR  ST. 

HARTFORD,  CT  06120 

METEOROLOGY  RESEARCH,  INC. 

464  W.  WOODBURY  RD. 

ALTADENA,  CA  91001 

CONTROL  DATA  CORP. 

METEOROLOGY  DEPT.  RSCH.  DIV. 
2800  E.  OLD  SHAKOPEE  RD. 

BOX  1249 

MINNEAPOLIS,  MN  55440 

OCEAN  DATA  SYSTEMS,  INC. 

2460  GARDEN  ROAD 
MONTEREY,  CA  93940 

DIRECTOR  OF  RESEARCH  (2) 
ADMINISTRATION 
CODE  012 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 

PROFESSOR  H.  FREPRICKSEN 
CHAIRMAN,  DEPT  OF  MATHEMATICS 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 

DR.  RICHARD  LAU 
OFFICE  OF  NAVAL  RESEARCH 
800  QUINCY  ST. 

ARLINGTON,  VA  22217 

PROFESSOR  G.M.  NIELSON 
DEPT.  OF  COMPUTER  SCIENCE 
ARIZONA  STATE  UNIVERSITY 
TEMPE,  AZ  85287 

DR.  H.  JEAN  THIEBAUX 
NATIONAL  METEOROLOGICAL  CENTER 
DEVELOPMENT  DIVISION 
WASHINGTON,  DC  20233 

PROFESSOR  PETER  ALFELD 
DEPARTMENT  OF  MATHEMATICS 
UNIVERSITY  OF  UTAH 
SALT  LAKE  CITY,  UT  84112 


MR.  ROSS  HOFFMAN 
ATMOSPHERIC  AND  ENVIRONMENTAL 
RESEARCH,  INC. 

840  MEMORIAL  DRIVE 
CAMBRIDGE,  MA  02139 

DR.  WILLIAM  A.  HOPPEL 
OFFICE  OF  NAVAL  RESEARCH 
800  QUINCY  ST. 

ARLINGTON,  VA  22217 

DR.  DAVID  F.  PARRISH 

NATIONAL  METEOROLOGICAL  CENTER 

NWS 

NOAA 

WASHINGTON,  DC  20233 

CHIEF  OF  NAVAL  RESEARCH  (2) 
LIBRARY  SERVICES,  CODE  784 
BALLSTON  TOWER  #1 
800  QUINCY  ST. 

ARLINGTON,  VA  22217 

COMMANDING  OFFICER 
NOARL 

NSTL,  MS  39529 

SUPERINTENDENT 
LIBRARY  REPORTS 
U.  S.  NAVAL  ACADEMY 
ANNAPOLIS,  HO  21402 

NAVAL  POSTGRADUATE  SCHOOL 
OCEANOGRAPHY  DEPT. 

MONTEREY,  CA  93943 

COMMANDER  (2) 

NAVAIRSYSCOM 

ATTN:  LIBRARY  (AIR-7226) 

WASHINGTON,  DC  20361 

COMMANDER 

NAVOCEANSYSCEN 

DR.  J.  RICHTER,  CODE  532 

SAN  DIETO,  CA  92152 

DR.  JOHN  HOVERMALE 
DIRECTOR,  ATMOSPHERIC  SCIENCES 
DIVISION 

NAVAL  OCEANOGRAPHIC  AND 
ATMOSPHERIC  RESEARCH  LABORATORY 
MONTEREY,  CA  93943 
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SUPERINTENDENT 
ATTN;  USAPA  (DEG) 

COLORADO  SPRINGS,  CO  80840 
AFGL/LY 

HANSCOM  APB,  MA  01731 

OFFICER  IN  CHARGE 
SERVICE  SCHOOL  COMMAND 
DET.  CHANUTE/STOP  62 
CHANUTE  AFB,  IL  61868 

DIRECTOR  ( 2 ) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA  22314 

DIRECTOR 

NATIONAL  METEOROLOGICAL  CENTER 
NWS,  NOAA 
WWB  W32,  RM  204 
WASHINGTON,  DC  20233 

DIRECTOR 

OFFICER  OF  PROGRAMS  RX3 
NOAA  RESEARCH  LAB 
BOULDER,  CO  80302 

DIRECTOR 

NATIONAL  SEVERE  STORMS  LAB 
1313  HALLEY  CIRCLE 
NORMAN,  OK  73069 

DIRECTOR 

TECHNIQUES  DEVELOPMENT  LAB 
GRAMAX  BLDG. 

8060  13TH  ST. 

SILVER  SPRING,  MD  20910 

EXECUTIVE  SECRETARY,  CAO 
SUBCOMMITTEE  ON  ATMOS.  SCI. 
NATIONAL  SCIENCE  FOUNDATION 
RM.  510,  1800  G.  STREET,  NW 
WASHINGTON,  DC  20550 

ATMOSPHERIC  SCIENCES  DEPT. 

UCLA 

405  HILGARD  AVE. 

LOS  ANGELES,  CA  90024 

CHAIRMAN,  METEOROLOGY  DEPT. 
CALIFORNIA  STATE  UNIVERSITY 
SAN  JOSE,  CA  95192 


NATIONAL  CENTER  FOR  ATMOS. 

RSCH.,  LIBRARY  ACQUISITIONS 
P.O.  BOX  3000 
BOULDER,  CO  80302 

CHAIRMAN,  METEOROLOGY  DEPT. 
PENNSYLVANIA  STATE  UNIVERSITY 
503  DEIKE  BLDG. 

UNIVERSITY  PARK,  PA  16802 

UNIVERSITY  OF  HAWAII 
METEOROLOGY  DEPT. 

2525  CORREA  ROAD 
HONOLULU,  HI  96822 

ATMOSPHERIC  SCIENCES  DEPT. 
OREGON  STATE  UNIVERSITY 
CORVALLIS,  OR  97331 

CHAIRMAN 

ATMOS.  SCIENCES  DEPT. 
UNIVERSITY  OF  VIRGINIA 
CHARLOTTESVILLE,  VA  22903 

CHAIRMAN,  METEOROLOGY  DEPT. 

UNIVERSITY  OF  UTAH 

SALT  LAKE  CITY,  UT  84112 

ATMOSPHERIC  SCI.  RSCH.  CENTER 
NEW  YORK  STATE  UNIVERSITY 
1400  WASHINGTON  AVE. 

ALBANY,  NY  12222 

METEOROLOGY  INTL. ,  INC. 

P.O.  BOX  22920 
CARMEL,  CA  93922 

SCIENCE  APPLICATIONS,  INC. 

205  MONTECITO  AVENUE 
MONTEREY,  CA  93940 

EUROPEAN  CENTRE  FOR  MEDIUM 
RANGE  WEATHER  FORECASTS 
SHINFIELD  PARK,  READING 
BERKSHIRE  RG2  9AX,  ENGLAND 

DEPT.  OF  MATHEMATICS 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 

PROFESSOR  R.  FRANKE,  MA/FE 
DEPT.  OF  MATHEMATICS  (10) 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 
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PROFESSOR  R.E.  BARNHILL 
DEPT.  OF  COMPUTER  SCIENCE 
ARIZONA  STATE  UNIVERSITY 
TEMPE,  AZ  85287 

LIBRARY,  FLEET  NUMERICAL  (2) 
OCEANOGRAPHY  CENTER 
MONTEREY,  CA  93943 

PROFESSOR  GRACE  WAHBA 
DEPARTMENT  OF  STATISTICS 
UNIVERSITY  OF  WISCONSIN 
MADISON,  WI  53705 

DR.  R.S.  SEAMAN 
AUSTRALIAN  NUMERICAL 
METEOROLOGY  RESEARCH  CENTRE 
P.O.  BOX  5089AA 
MELBOURNE,  VICTORIA, 

AUSTRALIA,  3001 

PROFESSOR  MARK  E.  HAWLEY 
DEPARTMENT  OF  ENVIRONMENTAL 
SCIENCES 

UNIVERSITY  OF  VIRGINIA 
CHARLOTTESVILLE,  VA  22903 

PROFESSOR  JAMES  J.  O'BRIEN 
DEPARTMENT  OF  METEOROLOGY 
THE  FLORIDA  STATE  UNIVERSITY 
TALLAHASSEE,  FLORIDA  32306 

NAVAL  OCEANOGRAPHIC  AND  (6) 

ATMOSPHERIC  RESEARCH  LABORATORY 
ATTN:  DR.  EDWARD  BARKER 

DR.  THOMAS  ROSMOND 
DR.  JAMES  GOERSS 
NANCY  BAKER 
PAT  PHOEBUS 
LARRY  PHEGLEY 
MONTEREY,  CA  93943 

DR.  LICIA  LENARDUZZI 
CNR- I AMI 

VIA  A.  M.  AMPERE,  56 
20131  MILANO 
ITALIA 


DR.  LEV  GANDIN 

NATIONAL  'meteorological  CENTER 
DEVELOPMENT  DIVISION 
WASHINGTON,  DC  20233 

DR.  ANDREW  LORENC 
BRITISH  MET.  OFFICE 
BRACKNELL  RB12  2SZ 
BERKS . ,  ENGLAND 

DEAN  G.  DUFFY 

NASA/GODDARD  SPACE  FLIGHT  CENTER 
LABORATORY  FOR  ATOMOSPHERES 
GREENBELT,  MD  20771 


COMMANDING  OFFICER  ( 2 ) 

NAVAL  OCEANOGRAPHIC  AND 
ATMOSPHERIC  RESEARCH  LABORATORY 
ATTN;  DR.  MOSELY 
NSTL,  MS  39529 


DR.  EMMETT  MADDRY 
K-43 

NAVAL  SURFACE  WARFARE  CENTER 
DAHLGREN,  VA  22448-5000 


DR.  LAWRENCE  C.  BREAKER 
NATIONAL  METEOROLOGICAL  CENTER 
MARINE  PREDICTION  BRANCH 
WASHINGTON,  DC  20233 
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